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SECTION  I 
INTRODUCTION 

Laser-Induced  damage  to  transparent  optical  materials  Is  a process  by  which 
optics  suffer  an  observable  change  upon  passage  of  a laser  beam.  Laser-Induced 
damage  can  be  divided  Into  four  temporal  regions:  damage  Involving  picosecond 
pulses,  pulses  from  a few  nanoseconds  to  a few  hundred  nanoseconds,  pulses 
from  a few  to  about  20  microseconds,  and  CW  laser  beams.  Over  the  longer 
periods,  damage  Is  characterized  by  the  cumulative  effect  of  absorption  of  the 
laser  beam  and  subsequent  distortions  or  heating  of  the  material  to  Its  melting 
point.  In  the  nanosecond  and  picosecond  time  regions  the  onset  of  damage  Is 
similar  to  the  exponential  growth  of  Ionization  and  resultant  plasma  formation 
associated  with  electrical  breakdown.  It  Is  In  the  few  tens  of  nanosecond 
pulse  region  with  which  this  report  Is  concerned. 

The  purpose  of  this  study  was  to  examine  the  material  properties  and  fabri- 
cation techniques  which  directly  bear  on  the  resistance  of  optical  materials 
to  laser-induced  damage.  In  particular  the  properties  and  manufacturing 
techniques  of  thin  films  applied  to  optical  substrates  were  examined.  A prime 
purpose  of  this  study  was  to  correlate  the  observed  behavior  with  material 
properties  such  as  refractive  Index,  condition  of  the  surface,  and  material 
structure. 

In  any  high-power  laser  system  employing  pointing  and  tracking  or  beam 
colli motion,  lenses,  windows,  and  Isolators  must  Intercept  the  beam  before  the 
output  can  be  utilized.  With  a nominal  4 percent  reflection  at  each  surface 
only  17  surfaces  are  required  to  reduce  the  available  power  by  50  percent.  Since 
optical  systems  can  easily  contain  this  many  surfaces  It  becomes  necessary  to 
place  on  the  bare  components  anti  reflective  thln-fllm  coatings.^  The  surviv- 
ability of  thln-fllm  coatings  to  laser  beams  then  Is  of  the  utmost  Importance 
In  the  reliability  and  efficiency  of  high-power  laser  systems.  By  Increasing 
the  threshold  to  laser-induced  damage  of  an  optical  surface  one  can  reduce  the 
size  of  that  optical  component  by  Increasing  the  power  density  which  It  can  with 
stand.  Such  size  reduction  and  resulting  weight  reduction  are  very  Important  to 
any  laser  devices  which  are  designed  for  airborne  use. 
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Because  it  is  presently  uncertain  what  the  final  design  requirement-  of  a 
laser  weapon  system  will  be  in  terms  of  output  wavelength,  the  resulting 
material > problem  may  not  yet  be  solved.  Those  materials  which  are  appropriate 
at  one  wavelength  may  be  entirely  unsuitable  at  half  or  twice  that  wavelength. 
Because  of  this  uncertainty  it  is  imperative  that  a fundamental  understanding 
be  gained  of  those  material  properties  and  surface  finishing  techniques  as  well 
as  thin-fllm  characteristics  which  affect  the  surface  threshold  to  laser-induced 
damage.  Due  to  surface  Imperfections  and  contamination  it  is  the  surface  of 
materials  which  generally  fall  first.  Surface  thresholds  are  generally  only 
0.3  to  0.5  of  the  bulk  damage  threshold.  Thus  not  only  material  properties  but 
the  mechanics  of  surface  finishing  determine  the  damage  threshold  of  optical 
elements  and  limit  the  power  density  which  can  be  accepted  by  the  element  with- 
out degradation.  Thus,  It  Is  frequently  necessary  to  expand  the  beam  and  optics 
to  lower  the  Incident  power  density  while  maintaining  the  total  transmitted 
power.  This  causes  a multitude  or'  problems.  Figuring  a 20- Inch  window  to 
X/10  tolerances  Is  more  expensive  and  more  difficult  to  do  than  a 10-Inch 
window.  Cooling  problems  climb  rapidly  with  diameter  and  the  resulting 
weight  penalties  and  Increased  volume  requirements  may  outstrip  the  ability 
of  a transport  vehicle  to  accommodate  It.  Thus,  It  Is  seen  that  the  phenomenon 
of  laser-induced  damage  to  optical  components  places  constraints  In  the  design 
of  high-power  laser  systems.  By  Increasing  the  resistance  to  laser-induced 
damage  of  the  materials  used  In  a high-power  laser  system  one  can  reduce  the 
system  In  weight  and  volume,  thereby  reducing  aerodynamic  Influences.  Con- 
comitant gains  In  operation  can  be  realized  since  slew  rates  can  be  made  greater 
for  smaller  systems  and  thus  pointing  and  tracking  can  be  Improved.  In  finding 
ways  to  Improve  the  damage  resistance  of  optical  materials  one  must  first  iso- 
late as  much  as  possible  those  factors  which  affect  the  damage  threshold.  With 
the  exception  of  laser-induced  damage  In  the  bulk  of  materials  caused  by  non- 
linear self  focusing  It  Is  the  surface  of  materials  which  falls  first.  An 
Important  question  to  be  considered  In  this  study  Is  which  factors  affect  the 
threshold  power  density  for  laser-induced  surface  damage.  Since  thin-fllm 
coatings  are  used  not  only  as  anti  reflective  devices  but  also  to  afford  a mea- 
sure of  environmental  protection  for  certain  materials.  It  Is  appropriate  to  ask 
why  films  fall  at  lower  Incident  laser  power  densities  than  do  bare  surfaces  of 
the  same  material.  In  seeking  an  answer  to  this  question,  one  Is  led  to  study 
those  factors,  both  Intrinsic  and  extrinsic  to  the  thin- film,  which  affect 
damage  thresholds.  These  considerations  place  a need  on  Identification  of  the 
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basic  mechanisms  involved  in  laser-induced  damage  so  that  functional  relation- 
ships between  material  properties  and  surface  treatment  and  damage  threshold 
levels  will  allow  selection  and  development  of  materials  which  will  meet  these 
engineering  requirements. 

A comprehensive  compendium  of  papers  dealing  with  laser-induced  damage  would 
probably  serve  no  useful  purpose.  The  subject  Is  so  Important  that  an  annual 
symposium  discussing  only  laser-induced  damage  and  attended  by  150  or  more 
scientists  to  hear  over  40  select  papers  is  now  In  Its  seventh  year.  During 
the  initial  stages  of  a new  science  or  the  Investigation  of  a new  process  It 
Is  difficult  for  the  experimentalist  to  identify  the  Important  parameters 
which  must  be  reported  to  completely  specify  the  experiment.  A great  number 
of  the  early  papers  on  laser-induced  damage  failed  to  report  such  factors  as 
pulse  widths,  focusing  conditions,  beam  quality,  characterization  of  meaningful 
material  properties  and  experimental  technique.  It  appears  now  that  each  of 
these  is  Important  If  one  were  to  attempt  to  reproduce  a given  laser  damage 
experiment  or  attempt  a correlation  of  published  data.  Omitting  these  reports, 
which  in  light  of  subsequent  developments  failed  to  adequately  specify  the 
experiment,  and  eliminating  those  which  deai  with  CM  systems,  picosecond  systems 
(except  as  relative  Indicators),  and  10.6  um  systems  leaves  primarily  those 
studies  conducted  and  reported  since  the  Inception  of  the  Symposium  on  Laser- 
Induced  Damage  to  Optical  Components  (refs.  1 through  7). 

Most  publications  on  laser-induced  damage  list  threshold  values  In  terms  of 
energy  density;  a few  list  damage  thresholds  In  terms  of  Incident  power  density; 
however,  an  elegant  and  simple  description  by  Crisp  et  el.  (ref.  8)  has  shown 
that  the  correct  quantity  in  specifying  the  damage  threshold  in  the  short  pulse 
region  is  the  electric  field  associated  with  the  laser  pulse.  This  proposal 
stems  from  the  well-known  observation  that  optical  flats  exhibited  lower  damage 
threshold  In  terms  of  energy  and  power  density  at  the  exit  face  than  at  the 
entrance  face  (ref.  9).  Many  complex  explanations  were  put  forth,  but  It  was 
Crisp  who  assumed  that  the  damage  threshold  depended  on  the  optical  electric 
field  and  that,  upon  applying  the  appropriate  Fresnel  equations,  the  effective 
electric  field  at  the  exit  face  of  an  optical  flat  was  greater  than  that  at  the 
entrance  for  a given  Incident  power  density  in  just  the  right  proportion  to 
account  for  the  apparently  lower  energy  threshold.  In  fact,  the  threshold 
electric  field  at  the  entrance  and  exit  were  the  same  (ref.  10).  As  proof,  an 
optical  flat  was  placed  In  a damaging  laser  beam  at  Brewster's  angle  so  that 
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reflection  losses  were  eliminated;  it  was  observed  that  the  entrance  and  exit 
faces  failed  at  the  same  level  of  laser  <-'-^nsity. 

In  the  work  of  DeShazer  et  al.,  and  Newman  and  DeShazer  (refs.  11,  and  12) 
further  support  Is  found  for  tying  the  electric  field  to  laser-induced  damage 
thresholds.  Films  of  Ti02  In  A/4,  A/2,  and  3A/4  thickness  were  deposited  on 
glass  substrates  and  subjected  to  damaging  laser  pulses.  The  threshold  fields 
were  5.13  x 107,  3.78  x 107  ana  3.69  x 107  V/m,  respectively.  These  values 
are  sensibly  constant  within  the  quoted  experimental  uncertainty  of  ±15  percent. 
These  values  are  apparently  related  more  closely  to  an  intrinsic  property  of 
T1 02  than  are  the  widely  variant  incident  threshold  power  densities  of  3.45, 

0.6,  and  1.8  GW/cm2.  This  divergence  in  threshold  power  densities  was  enhanced 
by  the  wide  variation  in  Index  between  the  film  (n  * 2.28)  and  glass  substrate 
(n  * 1.51).  The  relationship  between  refractive  Index,  laser  power  density,  and 
optical  electric  field  Is  contained  In  the  formulation  of  the  Poyntlng  Theorem 
taking  Fresnel  reflections  into  account.  The  Poyntlng  vector,  £(watts/m2),  Is 
related  to  the  electromagnetic  components,  t and  ft,  of  the  light  wave  by 

5 » t x fl  (1 ) 

The  constituent  relations  yield  the  ratio  between  electric  and  magnetic  field 
amplitude*  as 

Ve  E - V V H (2) 

where  e and  u are  the  dielectric  permittivity  and  permeability,  respectively. 

For  a linearly  polarized  plane  wave  In  a linear,  nonmagnetic  material  the 
magnitude  of  the  Poyntlng  vector  is  related  to  the  magnitude  of  the  optical 
electric  field  by 


The  identities  n2  ■ z/tQ  an.;  c2  * l/e0u0,  where  n and  c are  the  refractive 
index  of  the  medium  and  the  speed  of  light  In  fr«  ace,  respectively,  allows 
one  to  rewrite  equation  (3)  as 
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(4) 


Here  the  refractive  index  of  air,  nQ  * 1 , has  been  used  to  Indicate  the  magni- 
tude of  the  electric  field  incident  at  the  air-dielectric  interface.  The 
Fresnel  relations  are  now  used  to  calculate  the  optical  electric  field  at  the 
surface  of  a bare  substrate  and  at  the  surface  of  thin  films  of  A/4,  A/2,  and 
3A/4  thicknesses.  The  optical  electric  field  at  the  surface  of  an  uncoated 
material  with  Index  ng  Is 


E(V/m)  - -£4Vs(watts/m*)  (5) 

g 

For  films  of  A/4  and  3A/4  optical  thickness  with  Index  n^.  applied  to  a sub- 
strate of  index  ng,  the  optical  electric  field  at  the  film  surfaces  are 


E(V/m) 


'V 


nf 


19.4  Vs (watts 'm2) 


(6) 


For  a film  of  A/2  optical  thickness  the  field  is 

E(V/m)  - t-+”7T  19-4  Vs(watts/m2)  (7) 

i rig 

As  n^  approaches  ng  the  expressions  given  In  equations  (6)  and  (7)  coalesce  and 
if  the  field  is  indeed  the  threshold  quantity  one  would  expect  the  variation  In 
threshold  power  density  with  film  thickness  to  vanish.  Correspondingly,  we  note 
that  for  the  case  of  S102  (n  ■ 1.45)  applied  to  BSC-2  (n  » 1.51)  the  values  for 
threshold  field  versus  Incident  power  density  for  A/4  and  A/2  films  were  1.65, 
1.56  x 10*  V/m;  and  10.5,  9.9  GW/cm2,  respectively. 

Two  trends  apparent  from  a review  of  past  work  are  that,  In  general  break- 
down, threshold  varies  inversely  with  refractive  index.  From  Turners'  work 
(ref.  13),  In  which  he  tested  13  quarter  wave  films  at  0.69  un»,  it  was  seen  that 
damage  thresholds  were  lower  for  high  Index  materials  as  compared  to  low  Index 
materials.  This  fact  was  noted  by  Austin  et  al.  (ref.  14)  when  they  chose 
candidate  material r to  test  in  vapor-phase  mixture  films.  Their  work  also 
displayed  the  Inverse  dependence  of  damage  threshold  on  refractive  Index. 
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The  thin-fllm  studies  of  DeShazer  et  al.,  and  Newman  and  DeShazer  (refs.  11  and 
12)  displayed  a similar  trend  as  did  the  surface  studies  on  uncoated  optics  con- 
ducted by  Boling  et  al.  (ref.  15).  Because  surface  damage  studies  and  thin- 
fllm  studies  involve  many  more  variables  than  do  bulk  damage  studies  the  results 
of  Fradin  and  coworkers  (refs.  16  and  17)  were  taken  as  more  significant  of  the 
refractive  index  trend  than  were  the  studies  mentioned  above.  In  their  experi- 
ment they  damaged  11  alkali  halide  crystals  by  using  pulsed  lasers  with  10.6-wm, 
1.06-ym,  and  0.69-ym  wavelengths  and  compared  the  threshold  electric  field  with 
reported  dc  dielectric  breakdown  thresholds.  Not  only  was  a definite  trend 
between  the  damage  threshold  and  refractive  Index  established,  but  from  the 
close  correlation  between  the  dc  breakdown  fields  and  the  jamaging  optical 
fields  (being  about  1.5  times  the  dc  fields),  a confirmation  of  the  importance 
of  the  optical  electric  field  In  the  damage  process  was  evident. 

Although  not  quantified.  It  was  known  at  least  from  1970  (ref.  9)  that 
surface  roughness.  Impurities,  and  Imperfections  reduced  the  threshold  of  sur- 
faces below  that  of  the  bulk.  Although  the  attendent  problems  of  self-focusing 
and  Inaccurately  known  focal  volumes  In  the  bulk  of  materials  make  the  compari- 
son between  bulk  and  surface  threshold  tenuous.  It  was  generally  reported  that 
surface  thresholds  were  only  0.3  to  0.5  times  that  of  the  bulk  thresholds.  This 
problem  was  attacked  In  a paper  by  Bloembergen  (ref.  18)  In  which  he  used  static 
field  theory  to  calculate  the  field  enhancement  expected  In  certain  dielectric 
cavities  which  model  surface  Imperfections.  Electric  field  enhancement  by 
factors  of  up  to  the  square  of  the  refractive  Index  were  attained  In  cavities 
of  not  unreasonable  size  and  shape.  This  Infers  a maximum  field  enhancement  of 

-2.25  for  S10, , for  example.  One  main  conclusion  of  his  work  was  that  surface 
2 0 
features  with  characteristic  dimensions  less  than  about  100  A (10  nm)  would  be 

unimportant  In  laser-induced  damage.  Bloembergen  argued  that  diffusion  of 

electrons  out  of  a volume  with  a 100  X characteristic  dimension  would  be  so 

rapid  that  pulses  of  a few  nanoseconds  would  have  to  be  increased  in  power 

density  to  that  of  Intrinsic  threshold  levels  to  compensate  for  the  m1grc..,on. 

There  have  been  studies  to  confirm  Bloembergen's  basic  hypothesis  that  a surface 

' a 

with  a finish  which  has  no  Irregularities  larger  than  100  A will  have  the  same 
threshold  as  the  bulk.  There  has  also  been  conjecture  that  the  100  A criterion 
was  too  stringent.  Fradin  and  Bass  (ref.  19)  super  polished  fused  silica, 
sapphire,  and  BSC-2  glass  in  an  attempt  to  raise  the  surface  threshold  to  that 
of  the  bulk.  With  conventional  polish  the  ratio  of  surface  to  bulk  thresholds 
ranged  from  0.5  for  sapphire  to  0.77  for  FTC-2  and  fused  silica.  After 
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polishing  by  a technique  which  recirculates  the  polishing  slurry  (bowl  feed 
polishing),  they  were  able  to  raise  the  surface  threshold  to  that  of  the 
bulk  for  BSC-2  and  fused  silica.  The  sapphire,  notoriously  difficult  to  polish, 
was  not  improved  by  this  technique.  No  other  study  has  claimed  to  have 
observed  surface  threshold  equal  to  that  of  the  bulk. 

With  this  background  It  was  decided  to  investigate  the  dependence  of  the 
threshold  optical  electric  field  for  laser-induced  damage  on  such  material 
properties  as  refractive  Index,  number  density,  and  structure.  Realizing  the 
Importance  of  surface  finish,  very  careful  finishing  procedures  and  surface 
measurements  were  required  of  the  vendor.  For  the  thin-fllm  portion  of  the 
study  optically  half  wavelength  thick  films  (at  1.06  pm)  were  applied  to 
carefully  i-'lnlshed  flats  primarily  of  fused  silica.  The  substrates  were 
polished  as  smooth  as  possible  In  large  batches  to  minimize  variation.  Film 
deposition  was  carefully  controlled  as  were  a number  of  extrinsic  variables  such 
as  background  pressure  in  the  coating  chamber  to  obtain  films  as  free  from 
structural  Inhomogeneities  as  possible. 

This-  report  Is  divided  Into  five  sections.  Following  the  Introduction 
are  Experiment,  Theory,  Results,  and  Conclusions.  The  general  plan  of  attack 
was  to:  (1)  obtain  a maximally  characterized  laser  beam  which  exhibited  minimal 

spatial,  temporal,  and  transverse  mode  variation,  (2)  use  the  laser  to  damage 
a set  of  carefully  manufactured  coated  and  uncoated  optical  flats  to  obtain 
threshold  damage  fields  versus  refractive  index,  surface  structure,  film 
structure  to  include  Inhomogeneities,  film  deposition  techniques,  and  atomic 
number  density,  (3)  develop  theoretical  expressions  sufficient  to  explain  the 
observed  behavior. 


! 
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SECTION  II 
EXPERIMENT 

In  the  experimental  phase  of  this  study  various  optical  samples  with  and 

+3 

without  film  coatings  were  subjected  to  1.06-ym  laser  radiation  from  a Nd 
in  glass  laser  controlled  to  operate  In  the  TEM00  transverse  mode.  The  laser 
used  was  a Compagnie  General  Electric  Model  VD640  which  has  normal  operating 
characteristics  of  35  nsec  pulse  width  and  up  to  500  joules  output.  A greatly 
reduced  version  of  this  system  was  used  to  produce  the  high  beam  quality  re- 
quired with  attendant  characteristics  of  a 30  to  40  nsec  pulse  width  and  -1 
joule  output.  Of  the  seven  stages  In  the  total  system  only  the  oscillator  and 
first  amplifier  were  retained.  Each  consisted  of  16-mm  diameter  rods;  the 
oscillator  being  300  mm  long  and  the  amplifier  500  mm  long.  The  oscll’ator, 
Q-swItched  by  a Pockel's  cell,  was  constrained  to  the  TEMQo  transverse  mode 
using  an  Intracavity  aperture  of  1.9-mm  diameter  placed  between  the  oscillator 
rod  and  output  mirror.  (;or  the  piano-piano  cavity  employed  in  the  oscillator, 
it  Is  difficult  to  obtain  a reliable  calculation  for  the  size  aperture  necessary 
to  constrain  the  laser  to  the  TEMQ0  mode.  Because  of  this,  a number  of  aper- 
tures ranging  from  about  1 mm  In  diameter  to  about  3 mm  In  diameter  was  fabri- 
cated, and  the  proper  one  was  selected  on  the  basis  of  resulting  far-fleld 
Intensity  distributions.  This  method,  although  time  consuming.  Indicated  a 
proper  selection  of  1.9  ± 0.05  urn.  The  1.8-mro  aperture  was  no  Improvement 
in  the  expected  far-fleld  Intensity  distribution  but  provided  significantly  less 
than  the  nominal  80  mj  output  available  from  the  1.9-mm  aperture. 

The  oscillator  output  was  amplified  by  the  single  stage  amplifier  to  a total 
output  of  380-410  mj.  With  the  503.5-mm  focal  length  lens  used  in  this  experi- 
ment the  laser  radiation  was  focused  to  a 147-vjjn  spot  size  to  the  e"2  power 
points.  The  resulting  optical  electric  field  was  6.7  MV/cm  In  air,  which  Is 
sufficient  to  damage  the  surface  and  bulk  of  most  conventional  dielectrics  In  a 
30  to  40  nsec  pulse. 

The  remainder  of  the  equipment  was  used  for  diagnostics  and  parameterization. 

Figure  1 Is  a schematic  diagram  of  the  optical  test  conditions.  The  laser 
beam  was  directed  down  an  optical  bench  with  the  aid  of  two  45®-100  percent 
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figure  1.  Experimental  Arrangenent 
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reflectors.  These  reflectors  were  gimbal  mounted  to  facilitate  alignment  of  the 
system.  The  samples  were  held  In  a translatlng-rotatable  mount  approximately 
12  meters  from  the  beam  waist  of  the  oscillator. 

The  laser  was  mounted  on  12  cm  x 15  cm  granite  rails  to  maintain  mechanical 
stability.  Due  to  this  arrangement  the  alignment,  checked  twice  each  working 
day,  was  remarkably  stable.  To  eliminate  fluctuations  In  ambient  conditions, 
the  room  In  which  the  test  facility  was  located  had  Its  own  environmental  con- 
trol system.  A separate  heating/air  conditioning  system  maintained  the  tem- 
perature to  ±2.0  degrees  while  a humidity  control  system  kept  the  relative 
humidity  at  45  percent.  An  air  filtration  system  removed  airborne  particulate 
matter.  The  laser  was  fitted  with  a dry  nitrogen  purge  system  which  aided  In 
keeping  the  rod  faces  dust  free. 

The  focusing  and  beam  attenuation  systems  were  mounted  on  a 2-meter  lathe 
bed  optical  bench.  This  bench  was  rigidly  attached  to  a 4-1 nch  thick  honeycomb- 
backed  steel  plate  for  stability. 

The  focusing  lens  was  a C‘"Aefully  selected  biconvex  fused  silica  lens  of  the 
highest  quality.  A startling  revelation  In  the  course  of  this  study  was  that 
very  expensive,  supposedly  high-quality  lenses  exhibited  large  variations  from 
sphericity.  The  lens  chosen,  503.5  mm  focal  length,  was  tested  along  with  a 
couple  of  dozen  other  candidates  with  a Foucault  knife-edge  device.  The  chosen 
lens  alone  showed  excellent  sphericity. 

The  laser  was  operated  at  constant  input  energy  to  maintain  the  temporal 
and  spatial  quality  of  the  output  beam.  The  energy  delivered  on  target  was 
controlled  by  attenuating  the  main  beam  with  5 mm  x 5 mm  nonsaturable  Schott 
optical  filters.  These  filters  were  placed  In  the  beam  at  a 6°  angle  to  the 
normal  to  prevent  regenerative  feedback.  This  necessitated  using  the  filters 
In  sets  of  two  or  four  to  minimize  steering  of  the  laser  beam.  The  maximum 
transmission  needed  in  this  study  was  35  percent  and  the  minimum  was  4.5  x 1(TZ 
percent.  The  set  of  64  filters  provided  variation  between  these  ranges  In 
steps  separated  In  total  transmission  by  10  percent. 

The  3.8-cm  diameter  samples  were  held  in  a mount  with  6 degrees  of  freedom. 
The  samples  were  canted  5 degrees  from  the  axis  of  the  laser  beam  to  prevent 
feedback  Into  the  laser.  Even  with  this  precaution  when  hlgh-transmlsslon  pulses 
caused  dense,  bright  plasmas,  the  reflection  from  these  plasmas  was  enough  to 
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generate  post-lasing  In  the  conventional  mode,  and  a great  deal  of  care  was 
taken  to  identify  and  prevent  this  situation  from  arising. 

In  a past  damage  symposium  (ref.  20),  it  was  noted  that  the  plasma  associ- 
ated with  laser  damage  disrupted  the  surface  of  an  optical  sample  to  several 
diameters  of  the  damaging  beam.  And  Indeed,  one  needs  only  to  condense  water 
vapor  on  the  surface  of  a damaged  sample  to  clearly  and  Incontrovertlbly  Iden- 
tify each  damage  site  and  the  surrounding  area  which  was  disrupted  by  the  plasma 
(more  on  this  to  follow).  To  confine  this  disrupting  influence  to  the  vicinity 
of  each  separate  damage  site,  a face-plate  aperture  mask  was  fabricated  and 
Installed  on  each  sample  prior  to  testing.  The  mask  was  drilled  with  fifty- 
four  2.0-mm  diameter  holes  to  allow  passage  of  the  damaging  147-um  diameter 
laser  beam.  The  center  1.25  cm  was  drilled  to  accommodate  a metallic  coating 
used  In  the  characterization  of  each  surface.  Each  2.0-nm  aperture  was  sur- 
rounded by  an  annulus  of  circuit  board  layout  tape  which,  when  pressed  onto  the 
face  of  the  sample,  prevented  leakage  of  the  plasma  and  leakage  of  ultraviolet 
radiation  to  surrounding  damage  sites.  Examination  of  damaged  samples  confirmed 
the  confinement. 

A spot  size  of  approximately  150  um  in  diameter  was  chosen  to  obviate 
erroneous  results  due  to  the  spot  size  dependence  of  laser  Induced  damage 
observed  by  DeShazer  et  al.,  (ref.  11).  In  their  study  of  laser-induced  damage 
to  thin  films,  they  noted  that  increasing  the  spot  size  decreased  the  damage 
threshold  until  a spot  size  greater  than  about  150  um  was  reached.  Above 
-150  um  diameter  there  was  no  spot  size  dependence.  The  result  was  explained  on 
the  ba?4s  of  an  Increased  probability  for  hitting  a defect  when  using  a larger 
spot  size.  Diagnostics  equipment  was  mounted  rigidly  either  to  the  optical 
bench,  steel  plate,  or  granite  rails.  The  diagnostic  tools  were  of  two  cate- 
gories: active  and  parametric. 

1.  ACTIVE  DIAGNOSTICS 

The  energy  and  temporal  shape  of  the  laser  pulse  was  constantly  monitored 
by  two  bl planar  photodiodes  with  S-1  photocathodes  and  by  a calibrated  pyro- 
electric energy  meter.  The  temporal  characteristics  of  the  oscillator  pulse 
were  monitored  by  photodiode  P0-1  shown  In  figure  1.  The  output  of  this  photo- 
diode was  displayed  on  a Tektronix  Model  519  oscilloscope  and,  through  the  first 
3000  shots,  was  recorded  for  each  event  by  a Polaroid  scope  camera.  Due  to  the 
minute  variation  In  this  pulse  ehtpe,  during  the  last  2800  shots  only  every 
fifth  shot  was  so  recorded.  Figure  2 Is  a trace  typical  of  the  oscillator 
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Figure  2.  Typical  Oscillator  Pulse 

pulse.  The  amplifier  output  was  monitored  by  a second  photodiode  also  shown  In 
figure  1.  In  addition  to  the  temporal  response,  this  photodiode  contained  an 
Integrating  circuit  which  detected  pre-  and  post-lasing  and  gave  a measure  of 
the  total  output  energy.  Because  conventional  lasing  of  microsecond  duration 
prior  to  the  giant  pulse  (pre-lase)  or  following  the  giant  pulse  (post-lase)  Is 
not  detectable  on  the  fast-sweep  oscilloscope  used  to  record  the  glar.t  pulse. 

It  Is  necessary  to  monitor  these  conditions  separately.  To  monitor  both  the 
giant  pulse  (tens  of  nanoseconds)  and  the  conventional  lasing  (microseconds),  It 
Is  necessary  to  Incorporate  an  Integrating  network  with  an  R-C  time  constant  of 
about  200  usee.  On  a typical  sweep  rate  of  200  usec/cm,  a giant  pulse  appears 
as  a step  (see  figure  3)  while  the  conventional  mode  lasing  (usee  pulses  for 
about  0.5  millisecond)  builds  steadily  with  time.  Any  regenerative  feedback 
caused  by  reflections  from  optical  elements  downstream  from  the  final  amplifier 
can  be  amplified  In  the  active  medium.  Such  feedback!  when  coupled  Into  the 
oscillator,  can  damage  the  optical  elements  In  the  oscillator.  The  feedback  can 
only  be  anpllfled  when  the  gain  Is  sufficient  and  the  Pockel's  cell  Is  open. 

It  appears  as  post-lasing  on  the  Integrated  photodiode.  Constant  monitoring 
of  the  output  of  the  Integrating  photodiode  gave  notice  of  this  condition  at 
which  time  downstream  optical  elements  were  slightly  misaligned  to  prevent  sub- 
sequent feedback.  Because  the  step  function  Is  an  Integrated  record  of  the 
light  Intensity  In  the  giant  pulse,  It  also  served  as  a crude  energy  monitor. 

The  energy  Incident  on  target  was  monitored  by  a pyroelectric  detector 
placed  after  the  beam  attenuating  filters  and  before  the  focusing  lens.  It 
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Figure  3.  Integrated  Photodiode  Trace  Exhibiting  P^e-  and  Post-Lasing 

has  a dynamic  range  from  1 joule  full  scale  to  20  mlcrojojles  full  scale. 
Readings  as  low  as  S percent  of  full  scale  were  reliable  to  ±4  percent  of  full 
scale.  Through  beam  splitting,  attenuation  at  the  detector  head,  and  reflection 
losses,  the  device  recorded  3.51  percent  of  the  energy  on  target.  On  the  two 
samples  which  exhibited  the  Icwcst  threshold  for  damage,  It  was  necessary  to 
change  this  ratio  to  record  these  lower  energy  pulses.  The  detector,  a Laser 
Precision  model  Rk  d230»  was  compared  to  a separate  Rk  3230  and  they  were  found 
to  track  to  a precision  of  ±4  percent.  The  detector  was  calibrated  against  two 
separate  carbon  c*me  thermal  calorimeters  (TRS  Model  117)  which  have  an  accuracy 
of  ±10  percent  compared  to  t silver  sphere  calorimeter.  The  rest  of  the  active 
diagnostics  were  used  to  define  the  occurrence  of  damage. 
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For  the  purpose  of  this  study  we  Have  defined  damage  as  any  Irreversible 
change  Indicated  by  one  of  four  detection  techniques:  (1)  visible  plasma  for- 

mation, (2/  laser-induced  scattering,  (3)  phase-contrast  microscopy,  and 
(4)  breath- fogging. 

Referring  to  figure  1,  two  observers  were  positioned  behind  1.06-um  blocking 
filters  at  positions  A and  B.  They  actively  observed  each  event  for  the  appear- 
ance of  the  bright  plasma  which  usually  accompanies  the  occurrence  of  laser- 
induced  damage.  In  about  98  percent  of  the  events  this  observation  was  suffi- 
cient to  positively  Identify  damage.  For  the  two  low  Index  materials,  ZnS  and 
ZnSe.  there  were  several  events  that  (even  In  a darkened  room)  failed 
to  oroduce  a visible  plasma  during  a damaging  event.  In  these  cases  a general 
brightening  was  usually  observed  on  the  surface  by  one  or  both  observers.  Be- 
cause of  the  qualitative  nature  of  this  visual  observation  of  damage.  It  was 
reinforced  by  the  three  other  techniques,  although  In  the  great  majority  of 
tests  the  detection  of  a plasma  proved  effective  as  a damage  Indicator. 

The  appearance  of  observable  scatter  from  a He-Ne  laser  Incident  on  the 
damage  site  was  used  as  a secondary  Indicator.  The  forward  scatter  of  He-Ne 
No.  1 at  position  B was  used  as  well  as  the  scatter  of  He-Ne  No.  2 observed  at 
position  A.  Since  the  latter  was  focused  on  the  front  surface  of  the  sample, 
any  change  in  the  pattern  viewed  on  a screen  placed  behind  position  A afforded 
a method  for  distinguishing  between  front  surface  damage,  back  surface  damage, 
and  low-lying  (l.e.,  near  the  front  surface)  bulk  damage  sites.  The  diffraction 
pattern  from  front  surface  damage  was  quite  distinctive.  Although  this  tech- 
nique proved  useful  In  many  cases,  the  final  determination  of  surface  damage  was 
made  by  the  most  sensitive  method  described  later,  l.e.,  breath-fogging. 

A third  detection  scheme  Involved  an  adaptation  of  Zernlka  phase  contrast 
microscopy  (ref.  21).  Briefly  stated,  when  a coherent  beam  of  light  Is  re- 
flected from  a scattering  surface  the  specular  background  component  Is  brought 
to  focus  by  a lens.  This  light  Is  Insensitive  to  the  gross  nature  of  the  sur- 
face. That  light  which  Is  scattered  out  cf  the  background  contains  components 
of  higher  spatial  frequency  and  Is  not  generally  brought  to  a sharp  focus  by  the 
lens.  By  suitably  altering  the  phase  of  the  central  order,  It  Is  possible  to 
transform  the  phase  change  Introduced  by  a scattering  site  Into  an  Intensity 
variation  of  the  scattered  beam.  The  damage  sites  produced  by  near- threshold 
laser  pulses  are  not  discernible  In  ordinary  room  light.  They  absorbed  little 
or  no  light  and  as  such  are  largely  transparent,  When  light  passes  through  or 
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Is  reflected  from  one  of  these  small  damage  sltas,  the  predominant  effect  Is  the 
Intraduction  of  a spatially  varying  phase  shift  which*  since  ordinary  optics 
responds  to  changes  In  light  Intensity.  Is  not  directly  observable.  Zemlke 
protvsed  a phase  contrast  technique  which  not  only  Increases  the  contrast  and 
hence  the  visibility  of  such  phase  objects,  but  under  certain  circumstances  the 
observed  Intensity  varies  linearly  with  the  phase.  To  understand  this  process, 
assume  that  a damage  site  Introduces  a phase  shift  $(x,y)  into  coherent  light 
transmitted  through  It.  In  the  physical  layout  of  this  technique  the  coherent 
light  was  reflected  from  the  site,  but  this  merely  Introduces  a reflection 
coefficient  Into  what  follows  and  was  dropped  for  simplicity.  Following  Goodman 
(ref.  21)  the  light  transmitted  through  the  site  then  has  an  electric  field 
solitude  based  on  unit  Input  of 

u(x,y)  » exp[j$(x,y)]  (8) 

For  small  variations  In  phase  this  becomes 

u(x^)  s 1 + Jd(x,y)  (9) 

neglecting  terms  of  order  and  higher.  The  first  term  in  equation  (9)  is 
just  the  field  that  one  would  observe  If  the  damage  site  were  not  present,  while 
the  second  term  generates  weaker  diffracted  light  that  Is  scattered  out  from  the 
axis  of  the  system. 

The  Intensity  I(x,y)  then  Is 

Kx,y)  » | u(x,y)  |l  = | i + j*(x,y)  |*  s 1 (10) 

or  just  the  Input  Intensity.  Since  the  background,  corresponding  to  the  unper- 
turbed Input  wave.  Is  brought  to  a sharp  focus  while  the  diffracted  light  in- 
troduced by  the  phase  change  of  the  damage  site  and  containing  higher  spatial 
frequencies  Is  spread  away  from  the  focal  point.  It  Is  possible  to  adjust  the 
phase  of  the  central  order  so  that  it  interferes  with  the  diffracted  light. 

By  Introducing  • phase  shift  of  e/2  radians  into  the  focused  background, 
one  obtelns 

I(x^)  « | exp( j*/2)  * to  I*  « | jO+q)  i*  s?  1 ♦ to  (11) 

21 


AFWt-TR-76-61 


Thus  the  Intensity  has  become  linearly  related  to  the  phase  shift  The  re* 

qulred  phase  shift  was  Introduced  by  depositing  a X/4  optically  thick  transpar- 
ent dielectric  dot  on  a glass  substrate  so  that  the  dot  just  intercepted  the 
focused  background  radiation.  It  was  hoped  that  this  technique  could  be  used 
to  provide  a quantitative  measure  of  the  radial  extent  and  depths  of  the  dis- 
ruptions at  each  site  to  be  able  to  compare  the  amount  of  material  removed  with 
the  Incident  energy  density.  To  this  end,  a comparison  was  made  between  the 
Information  Inferred  from  the  phase  contrast  method  and  physical  measurements  of 
selected  craters  using  a Talystep  weighted  stylus  surface  profiler.  A positive 
correlation  was  obtained*  which  did  not  lend  Itself  readily  to  quantitative 
analysis.  Figure  4 Is  one  scan  done  by  the  Talystep  on  a laser-produced  crater 
In  SF-12  glass.  The  different  materials  exhibited  crater  profiles  more  or  less 
characteristic  of  the  material.  Note  the  difference  In  figure  5*  which  Is  a 
Talystep  done  on  a crater  In  fused  silica.  Here  the  sides  do  not  drop  abruptly 
from  the  average  surface  but  gently  slope  for  several  radii  of  the  central 
crater.  Finally  figure  6 depicts  the  ring  structure  typical  of  damage  In  ED-2 
and  BK-7  glasses.  These  craters  were  produced  by  laser  pulses  with  power  densi- 
ties generally  well  above  threshold  but  were  Impossible  to  locate  without  the 
final  damage  detection  scheme. 

The  final  technique,  which  proved  the  most  sensitive,  was  used  only  after 
the  sample  was  removed  from  the  holder.  It  has  been  noted  (ref.  20)  that  con- 
densing water  vapor  on  the  damaged  surface  of  an  optical  sample  clearly  In- 
dicated the  extent  of  the  damage  site.  It  has  also  been  noted  (ref.  20)  and 
mentioned  above  that  the  plasma  associated  with  laser-induced  damage  alters  the 
surface  surrounding  the  damage  site.  During  the  course  of  this  study  these  two 
observations  were  tied  together.  Upon  receipt  of  samples  which  had  been  Ion 
polished,  a cleaning  procedure  was  Instigated  which  Included  as  one  step 
rinsing  In  water.  On  the  Ion-polished  samples  the  water  would  not  wet  the  sur- 
face but  beaded  up  and  was  removed  readily  by  agitation.  The  phenomenon,  what- 
ever it  Is,  must  account  for  the  efficacy  of  the  breath-fogging  technique. 
Plasmas  associated  with  damaging  laser  Irradiation,  In  effect,  cause  Ion  bom- 
bardment of  the  surface  In  the  vicinity  of  the  Interaction  zone.  When  one 
blows  his  breath  across  the  saqple,  the  first  places  that  the  condensed  water 
vapor  escapes  from  are  the  plasma-polished  damage  sites.  When  viewed  In 
oblique  white  11$it,  the  damage  sites  clearly  show  up.  Whenever  there  was  a 
doubt  about  whether  a site  had  sustained  damage,  the  fogging  technique  was 
employed.  Care  must  be  taken  when  employing  this  technique  to  ensure  that  no 
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Figure  4.  Talystep  of  Laser-Produced  Crater  In  SF-12  Glass 


Figure  5.  Talystep  of  Laser-Produced  Crater  In  Fused  Silica 


Figure  6.  Talystep  of  Laser-Produced  Crater  In  BK-7  Glass 
Showing  Ring  Structure 
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damage  testing  Is  attempted  subsequent  to  fogging  the  sample.  The  surface 
threshold  is  markedly  reduced  after  the  surface  has  been  fogged.  It  proved  to 
be  the  most  sensitive  of  the  techniques.  The  alteration  of  the  surface  by  the 
plasma  appears  to  be  permanent  sine,  samples  damaged  more  than  a year  ago  still 
reveal  their  sites  by  breath-fogging. 

2.  PARAMETERIZATION 

It  has  been  the  general  opinion  in  recent  years  that  a TEMQ0  laser  was 
essential  in  the  study  of  laser-induced  damage.  It  was  argued  (ref.  22)  that 
filamentary  hot  spots  in  the  beam  contained  a fraction  of  the  total  energy 
sufficient  to  cause  damage  yet  was  contained  in  such  a small  area  that  actual 
energy  density  on  target  due  to  this  filament  could  be  an  order  of  magnitude 
greater  than  the  average  energy  density.  Recently,  however,  argument  has  been 
made  (ref.  23)  that  a TEMQ0  beam  Is  not  necessary  because  filamentary  structure 
In  a laser  beam  is  focused  to  a spot  size  too  small  to  contribute  to  damage. 
Real-world  lasers  are  not  always  ooerated  in  the  TEM_  mode,  and  thus  a certain 
justification  for  using  multimode  lasers  in  damage  studies  is  made.  However, 
to  understand  the  processes  and  assign  meaningful  numbers  such  as  threshold 
optical  electric  fields  to  damage  tests,  it  is  necessary  to  employ  a well- 
defined  beam  with  constant  spatial  pattern.  To  this  end,  a TEMQ0  mode  laser 
was  used  in  this  study.  To  check  on  the  TEMQ()  nature  of  the  beam  a near-field 
scan  at  10  m from  the  oscillator  waist  and  a far-field  scan  at  the  focal  plan 
of  both  a 1.9-m  lens  and  the  503.5-mm  lens  used  in  the  surface  damage  experi- 
ment was  performed.  Since  the  fundamental  mode  Is  a Gaussian  plane  wave  In  the 
far  field  (i.e.,  at  a position  where  the  Fraunhofer  approximation  is  valid), 
the  Intensity  profile  will  be  Gaussian  in  spatial  extent.  A measure  of  this 
spatial  profile  was  obtained  In  two  ways.  In  the  first  arrangement  shown  In 
figure  7,  a photodiode  was  used  as  a receiver  for  the  total  laser  pulse.  A 
130-um  diameter  pin  hole  supplied  by  A.  Owyoung  (ref.  24)  was  used  to  sample 
the  focused  beam  across  Its  approximately  4-mm  diameter.  The  temporal  shape  of 
the  sampled  pulse  was  compared  to  that  of  the  full  beam  as  revealed  by  photo- 
diode number  1 and  no  variations  were  found.  The  deflection  of  the  oscilloscope 
trace  from  the  photodiode  used  to  observe  the  sampled  beam  was  ratloed  to  the 
deflection  of  the  oscilloscope  trace  from  the  photodiode  used  to  observe  the 
whole  beam.  The  variation  of  this  ratio  as  the  pinhole  was  stepped  across  the 
output  beam  gave  a measure  of  the  radial  Intensity  distribution  of  the  beam. 

Beam  center  was  located  u.id  both  a vertical  and  horizontal  scan  were  made.  By 
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placing  a lens  a distance  of  one  focal  length  In  front  of  the  pinhole,  the 
half  power  diameter  was  found  and  related  to  the  full  angle  beam  divergence  by 
d ■ fe,  where  f Is  the  focal  length  of  the  lens  used  (ref.  25).  For  a trans- 
mission of  99  percent  of  a Gaussian  beam  through  an  aperture,  one  requires 
d ■ 3w  where  w Is  the  Gaussian  spot  size  and  d Is  the  diameter  of  the  aperture. 
The  aperture  used  was  1.9  mm,  thus  the  full -angle  beam  divergence  becomes 
1.07  mrad.  The  measured  beam  divergence  of  1.34  to  1.58  mrad  agrees  within 
experimental  accuracy  to  the  expected  value  of  1.38  times  diffraction  limited 
beam  divergence.  This  fact,  coupled  with  the  Gaussian  fit  to  the  Intensity  of 
13  to  17  percent  supports  the  claim  of  TEMQ0  mode  operation.  The  procedure 
was  repeated  using  the  503.3-mm  lens  and  replacing  the  receiver  photodiode 
with  the  Rk  3230  energy  detector.  The  results  In  figure  8 are  plotted  against 
a Gaussian  for  comparison.  For  the  1.9-m  lens  the  fit  was  ±13.6  percent  while 
for  the  503.5  mm  lens  the  fit  was  ±17  percent.  In  the  latter,  the  beam  diver- 
gence was  calculated  to  be  1.09  mrad. 

It  was  necessary  to  know  the  focused  spot  diameter  of  the  laser  beam  since 
the  surface  of  the  optical  sample  was  placed  not  at  the  focal  plane  of  the 
503.5  mm  but  at  the  focus  of  530.3  mm  from  the  midplane  of  the  lens.  Since  a 
Gaussian  beam  remains  a Gaussian  under  propagation.  It  was  necessary  to  measure 
the  Gaussian  spot  size  at  the  focus.  The  method,  Illustrated  in  figure  9,  con- 
sists of  using  specially  treated  dark  paper*  which  bleaches  upor;  receiving  a 
specific  Incident  energy  density  threshold  of  * 2 j/cm2.  The  team  was  succes- 
sively attenuated  as  the  laser  was  fired  on  separate  sites.  For  a Gaussian  beam 
the  Intensity  varies  as  I (r)  as  I0e"r  /w*,  thus  a plot  of  (I / LQ)  versus  r2  gives 
a relationship  which  readily  Identifies  the  spot  size,  w.  By  this  method  It  was 
found  that  the  spot  size  to  the  e~2  power  point  was  147  pm.  A number  of  these 
determinations  throughout  the  course  of  the  study  agreed  to  better  than  10 
percent.  j 

3.  SURFACE' CHARACTERIZATION 

One  of  the  prime  reasons  for  conducting  this  study  was  to  determine  those 
characteristics  of  surface  finish  and  film  deposition  which  materially  affect 
the  laser- Induced  surface  threshold.  It  was  Imperative,  therefore,  to  fully 
characterize  the  surfaces  of  both  the  uncoated  samples  and  the  thln-fllm 
coatings. 
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The  surfaces  were  measured  for  roughness  by  three  methods:  fringes  of  equal 

chromatic  order  (FECO),  total  Integrated  scatter  (TIS),  and  Talystep  depth 
profiling.  FECO  was  the  prime  Indicator  of  surface  roughness  while  the  total 
Integrated  scatter  and  Talystep  were  used  only  as  auxiliary  measurements.  For 
a flat  surface  the  roughness  Is  defined  as  the  root  mean  square  deviation  of 
surface  heights  from  the  average  surface.  Thus  on  a surface  with  a specified 
rms  roughness,  surface  features  as  large  as  twice  that  In  depth  are  not  un- 
common. 

FECO  interferometry  Is  a technique  shown  schematically  In  figure  10.  Both 
a reference  flat  and  the  surface  to  be  measured  are  partially  silvered.  The 
sample  and  flat  are  brought  to  within  a wavelength  or  so  of  each  other  and 
parallel.  A portion  of  the  two  plates  is  Imaged  by  lens  Lz  onto  the  slit  of  a 
spectroscope.  As  the  distance  between  the  flat  and  sample  varies  due  to  surface 
roughness,  the  spectroscope  passes  only  those  wavelengths  which  satisfy  2d  * mX 
(m  an  Integer).  Since  m is  constant  along  any  one  fringe,  d/x  Is  constant.  The 
fringes,  referred  to  as  fringes  of  equal  chromatic  order,  reveal  the  topographi- 
cal features  of  the  surface  of  the  sample  Imaged  on  the  slit,  of  the  spectro- 
scope (refs.  26  and  27). 

Total  Integrated  scatter,  from  a discussion  by  Porteus  (ref.  28),  is  a 
technique  which  relates  the  scatter  from  the  reflection  of  a laser  beam  incident 
on  a surface  to  the  surface  roughness.  The  sample  is  normally  located  tangent 
to  one  port  of  a hollow,  diffuse  reflecting  sphere.  Laser  light  enters  another 
port  and  the  nonspectral  reflected  light  enters  a detector  after  being  totally 
Integrated  by  the  sphere  (figure  11).  Surface  roughness  Is  given  by 

o - ^ y-ln(l-TIS)  02) 

where  X Is  the  wavelength  of  the  test  laser  and  TIS  Is  the  fraction  of  the  test 
beam  scattered  outside  the  prime  spectral  reflection.  The  technique  also 
affords  a method  for  determining  scatter  of  imperfections  in  thln-fllm  coatings. 

The  Talystep  device  uses  a weighted  stylus  and  electronic  amplification  of 

Its  vertical  motion  to  topographically  scan  the  surface  of  the  sample.  For 

hard  films  it  can  also  be  used  to  mechanically  measure  the  physical  thickness 

of  the  film.  It  was  also  used  In  this  study  to  measure  a few  of  the  damage 

craters  (figures  4,  5,  and  6).  The  final  method  was  reliable  for  surfaces  with 

0 

values  of  a > „50  - 75  A rms,  but  the  diameter  of  the  tip  of  the  stylus  precluded 
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MIRROR  SPUTTER 


Figure  ' , Schematic  Representation  of  the  Total  Integrated  Scatter  (TIS) 
Method  for  Determining  Surface  Roughness 

the  mr"  ament  of  very  small  features  and  gave  uniformly  too  low  readings  on 

the  smoi  *r  samples.  The  smallest  stylus  radius  was  1.0  ym.  To  obtain  an  rms 
surface  ighness  by  Talystep,  a one-dimensional  scan  of  2 mm  in  length  Is 
obtained  at  several  representative  locations.  For  each  scan  of  the  topography 
of  the  surf?  e the  surface  height  Is  measured  from  a baseline  at  regular  Inter- 
vals to  det  mine  an  average  surface  height.  The  deviation  of  the  surface 
height  from  the  average  Is  then  measured  at  regular  Intervals  and  the  root  mean 
square  deviation  Is  calculated.  It  Is  this  deviation  which  Is  termed  surface 
roughness. 

From  the  results  of  surface  roughness  measurements  performed  by  the  vendor, 
the  author,  and  the  Optical  Evaluation  Facility  at  the  Naval  Weapons  Center, 
China  Lake,  California,  the  striking  conclusion  was  the  lack  of  correlation 
between  the  measurements.  For  example,  sanple  800-7169-080  was  supplied  by  the 
vendor  with  a roughness  of  42  l rms  as  determined  by  the  most  accurate  tech- 
nique, FECO.  Accompanying  that  reading  was  the  total  Integrated  scatter  of 
0,11  percent,  whleh  Indicates  an  equivalent  roughness  of  16.7  l rms.  The  same 
sample  was  measured  at  China  take  by  TIS  at  16. 5 I and  by  the  author  by 
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Talystep  at  27  A.  It  would  appear  that  the  state  of  the  art  In  surface  measure- 
ments still  allows  for  a great  deal  of  uncertainty  In  surface  roughness  deter- 
minations. Industry  sources  (ref.  29)  indicate  that  the  FECO  measurement  Is 
the  most  reliable,  although  very  time  consuming  and  difficult  to  perform  and 
analyze.  For  this  reason  the  FECO  measurement  was  used  In  calculations  in- 
volving surface  roughness. 

The  refractive  Indexes  of  the  bare  samples  and  films  were  measured  by 
elllpsometry.  Number  densities  and  nearest  neighbor  separations  were  calculated 
or  obtained  from  standard  references. 

Of  major  concern  In  the  study  of  the  thin  films  was  the  homogeneity  of  the 
film  throughout  Its  depth.  The  spectral  reflectivity  is  a reliable  Indicator 
of  this  parameter.  A Carey  14  spectrophotometer  was  used  In  the  double  pass 
reflection  mode  to  obtain  the  spectral  reflectivity  of  the  films.  Figure  12 
Is  the  spectral  response  of  two  A/2  ThF4  and  Zr02  films.  These  two  films  are 
examples  of  the  excellence  of  a film  coater's  art  In  being  able  to  produce 
homogeneous  coatings.  Since  both  films  were  specified  as  half  wave  at  1.06  urn, 
a perfect  film  would  have  the  reflectance  of  the  substrate  at  1.06  urn.  Since 
ThF4  at  n * 1.49,  matches  Index  more  closely  with  the  fused  silica  substrate 
index  at  1.449  than  does  the  ZrOa  at  n « 2.0,  the  reflectance  of  the  ThF^  film 
versus  wavelength  Is  not  very  different  from  the  substrate.  In  fact,  the  maxi- 
mum reflectance  of  the  ThF,,  film  on  S102  at  2.12  urn,  where  the  film  Is  a quarter 
wavelength  thick.  Is  only  4.15  percent  compared  to  3.36  percent  for  the  sub- 
strate, while  for  the  ZrOz  film  the  maxlmun  Is  22.3  percent  at  2.12  pm.  Con- 
versely figure  13  reveals  the  spectral  response  of  L1F  and  MgFa.  These  films 
display  a markedly  Inhomogeneous  character  with  an  index  apparently  Increasing 
outward  from  the  substrate.  Since  they  both  have  refractive  Indexes  less  than 
that  of  the  substrate,  their  reflectance  should  Increase  toward  and  match  that 
of  the  substrate  at  1.06  urn.  These  two  films  along  with  BaF2  exhibited  the 
most  Inhomogeneity  In  this  respect.  The  author  will  call  attention  to  this 
fact  during  the  discussion  of  the  thin- film  results. 

To  appreciate  the  magnitude  and  nature  of  the  Inhomogeneity,  a computer 
program  designed  by  Capt  John  Loomis  (ref.  30)  was  used  In  an  attempt  to  model 
the  observed  behavior  of  these  films.  The  program  was  designed  to  accept  up  to 
50  film  layers  of  varying  thickness,  absorptance,  and  refractive  Index  and  to 
compute  the  angular  and  wavelength  dependence  of  the  reflectance.  The  param- 
eters used  Included  the  total  physical  thickness  of  the  film  as  measured  by 
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Talystep,  the  refractive  Index  of  the  substrate,  and  the  measured  refractive 
Index  at  the  surface  of  the  film.  The  Inhomogeneous  film  with  varying  Index 
was  discretized  Into  a stack  of  from  four  to  eight  films  with  refractive  Index 
and  thickness  the  variable  parameters.  It  was  found  that  a variation  In  the 
refractive  Index  of  about  7 percent  (1.28  at  the  substrate  to  1.38  at  the  film 
surface)  produced  a spectral  response  quite  suggestive  of  the  measured  spectral 
response.  One  result  has  been  reproduced  as  figure  14  which  should  be  compared 
to  figure  13. 

If  one  assumes  that  fused  silica  Is  Invariant,  Independent  of  the  source, 
a comparison  of  the  absorption  and  transmission  spectra  for  various  fused 
silica  blanks  will  quickly  dispel  that  notion.  While  the  major  portion  of  this 
effort  used  fused  silica  supplied  by  the  Amersll  Corporation  under  the  trade 
name  Optosll  I,  an  ancillary  study  Involved  fused  silica  substrates  supplied  by 
a variety  of  manufacturers.  During  the  coating  process  absorption  spectra  are 
routinely  taken.  From  these  spectra  varying  amounts  of  hydroxyl  Ion  was  Identi- 
fied in  the  different  fused  silica  types.  This  observation  formed  the  basis 
for  an  unexpected  conclusion  which  at  first  seemed  to  be  anomalous  behavior. 

This  result  will  be  discussed  In  detail  later.  This  last  measurement  was  men- 
tioned to  Indicate  that  in  thln-fllm  coatings  It  is  not  at  all  entirely  obvious 
what  deposition  conditions,  techniques,  and  properties  of  films  lead  to  high 
damage  threshold  films. 

4.  SAMPLE  PREPARATION 

The  test  samples  for  this  study  required  surface  finishes  of  two  types: 

(1)  bare  surfaces  and  substrate  surfaces  for  film  depositions  with  the  lowest 
practical  surface  roughness  and  (2)  bare  surfaces  and  substrate  surfaces  for 
film  depositions  with  a controlled  range  of  surface  roughnesses.  The  latter 
category  was  designed  to  test  the  Influence  of  surface  roughness  on  the  laser- 
induced  damage  threshold  and  whether  thln-fllm  thresholds  are  affected  by  the 
surface  roughness  of  the  substrate.  Both  the  smooth  and  rough  surfaces  were 
created  by  variations  of  the  process  of  controlled  grinding. 

The  effect  of  the  grinding  process  In  the  preparation  of  optical  surfaces 
extends  below  the  surface.  A layer  of  subsurface  disorder  In  the  form  of  micro- 
cracks and  fissures  develops  when  an  optical  surface  Is  ground.  A coamonly 
accepted  rule  In  the  optical  Industry  (ref.  31)  Is  that  subsurface  disorder 
extends  to  a depth  roughly  three  times  the  diameter  of  the  grit  used  In  the 
grinding  process.  Controlled  grinding  Is  a technique  for  minimizing  the 
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subsurface  disorder.  The  controlled  grinding  process  used  to  obtain  the  "best 
finish"  surface  proceeds  as  follows:  A 30-ym  grit  Is  used  to  rough  In  the 
surface  finish.  Since  this  process  Introduces  a subsurface  disorder  to  a depth 
of  -90  urn  the  next  stage,  employing  20-um  grit,  proceeds  until  -100  urn  of  the 
material  Is  removed,  thus  eliminating  the  disorderly  effect  of  the  30-pm  grit. 
This  process  continues  through  12- , 5-,  and  3-um  diameter  alumina  grit,  each 
stage  removing  material  to  a depth  about  three  times  the  diameter  of  the  pre- 
vious grit.  The  surface  Is  then  polished,  which  removes  about  10  pm  of  mate- 
rial and  gives  the  final  finish.  The  final  roughness  values  varied  from  about 

O O 

10  A rms  for  fused  silica  to  about  40  to  50  A rms  for  sapphire  and  ZnS. 

The  roughness  samples  were  created  by  Interrupting  the  polishing  at  selected 
total  polish  times.  Samples  In  the  300  A(rms)  range  were  polished  for  about 
30  minutes.  Forty  minutes  of  polishing  time  gave  a 150  A(rms)  surface,  while 
10  to  20  hours  gave  40  to  50  A(rms),  and  best  finish  (less  than  20  A(rms)) 
required  40  hours  polishing  time. 

The  dielectric  films  were  deposited  by  one  of  three  frequently  employed 
techniques:  rf  sputtering,  electron  beam  heating,  and  thermal -evaporation.  A 
set  of  MgF2  samples  were  prepared  by  each  of  the  three  deposition  techniques. 

The  MgF2  films  applied  to  substrates  on  which  the  substrate  roughness  was 
varied  were  deposited  by  thermal  evaporation.  The  remainder  of  the  films  were 
applied  by  electron  beam  heating. 

Each  dielectric  film  was  specified  to  be  a half  wave  In  optical  thickness  at 
the  1.06-um  wavelength  used  In  the  study.  tThls  parameter  was  monitored  by 
observing  the  reflectance  of  the  film-substrate  system  as  the  film  was  being 
applied.  For  films  whose  Index  Is  greater  than  the  substrate  Index,  a minimum 
In  reflectance  Is  obtained  when  a half-wave  optical  thickness  Is  reached.  For 
homogeneous  films  this  minimum  corresponds  to  the  reflectance  of  the  bare  sub- 
strate. The  half-wave  thickness  was  chosen  because  this  assured  that  the 
electric  field  at  the  substrate- film  interface  was  the  same  as  the  optical 
electric  field  at  the  film-air  Interface.  It  was  also  the  thickness  which  gave 
the  least  variation  in  the  electric  field  of  a standing  wave  pattern  for  slight 
variations  In  the  thickness.  For  films  whose  Index  is  less  than  the  Index  of 
the  substrate,  a reflectance  maximum  Is  obtained  when  a half-wave  optical 
thickness  is  reached.  For  films  whose  Index  was  the  same  as  or  very  near  the 
Index  of  the  substrate.  It  Is  necessary  to  use  a separate  witness  plate  with  a 
sufficiently  different  refractive  Index  to  accurately  monitor  the  thickness. 
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The  materials  covered  In  figures  12  and  13  make  it  apparent  that  the  most  sensi 
tive  variation  of  reflectance  about  the  extremum  was  provided  by  the  Zr02 
(n  ■ 2.0)  on  fused  silica  (n  ■ 1.449).  Thus  to  obtain  the  most  sensitivity  in 
the  thickness  determination  requires  the  use  of  a witness  plate  whose  index 
varies  from  that  of  the  film  as  much  as  possible. 


s 
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SECTION  III 
THEORY 

Pulsed  laser-induced  damage  in  solid  dielectrics  is  intimately  linked  with 
the  study  of  the  dc  breakdown  of  solid  dielectrics.  This  assertion  follows 
from  the  electric  field  view  of  laser-induced  damage  which  now  has  wide  accep- 
tance for  the  reasons  noted  In  the  Introduction  of  this  report.  It  Is  not 
reasonable  to  claim  that  there  is  a one-to-one  correspondence  between  the  two 
processes.  Laser-induced  damage  does  not  Involve  electrodes  with  their 
associated  Influences  nor  Is  the  damage  region,  therefore  the  total  energy 
drop,  as  localized  In  dc  breakdown  studies  as  It  is  in  laser  damage  situations. 
It  is  appropriate  to  argue,  however,  that  certain  of  the  electric-field  related 
processes  are  similar  between  the  two  phenomena.  For  this  reason  a brief 
background  of  the  current  theories  of  dc  breakdown  of  solid  dielectrics  will 
be  given,  followed  by  a review  of  current  thought  on  laser-induced  damage  to 
solid  dielectrics,  and  finally  a proposed  model  and  theory. 

1.  DC  BREAKDOWN  OF  SOLID  DIELECTRICS 

A.  von  Hippel  (ref.  32)  In  an  attempt  to  reconcile  the  accepted  theory  of 
dc  gas  breakdown  to  solids  proposed  that  an  energy  balance  must  be  maintained 
for  any  electron  In  the  conduction  band  of  a solid.  A sufficient  alteration 
from  this  equilibrium  energy  balance  (e.g.,  by  applying  a large  enough  electric 
field)  leads  to  an  Increase  In  the  energy  of  every  free  electron  and  subsequent 
dielectric  breakdown.  Under  the  Influence  of  an  electric  field,  E,  the  z com- 
ponent of  momentum  will  obey 


W-)  E ' -«  (’3, 

However,  the  motion  of  the  electron  is  retarded  due  to  electron- phonon  colli- 
sions 
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where  t(E),  the  relaxation  time,  Is  assumed  to  be  a function  of  the  electron 
energy.  In  a steady  state  when 


(*) 


• 0 


(15) 


the  electron  drifts  In  the  field  direction  with  velocity 


m \m  / 


(E)E  - -u(E)E 


(16) 


where  m Is  the  electron  mass  and  u(E)  Is  the  mobility  of  the  electron.  This 
drift  velocity  Is  superimposed  cn  the  random  thermal  motion.  This  will  be  seen 
to  be  of  importance  later  when  we  consider  the  calculations  of  Seitz  (ref.  33). 
To  determine  the  total  energy  of  the  electron,  one  must  balance  the  gain.  A, 
experienced  from  the  applied  field  against  the  losses  to  the  lattice,  B. 

Energy  gain  Is  given  by 

q^EME) 

■ 07) 


A(E,E.T0)  - qeu(E)E2  - 


where  TQ  Is  the  lattice  temperature.  Electron- phonon  Interactions  Involve  the 
absorption  or  emission  of  a phonon  of  energy  hw.  The  ratio  of  the  probability 
of  emission  versus  absorption  is  given  in  elementary  quantum  theory  as 
(Nw+1)/Nw  where  Is  the  average  number  of  lattice  quanta  with  energy  hw, 
that  Is 


0) 


1 + exp 


(18) 


The  rate  at  which  energy  Is  lost,  B(E,T0),  Is  obtained  by  observing  that  a net 
loss  occurs  In  only  one  of  (2Nu  + 1)  collisions  so  that 


B(E,T0)  , - 


wfanTT 

* UJ 


(’  ■'r) 


(19) 
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where  tt(E)  Is  the  mean  time  between  collisions  and  Is  equal  to  t(E)  If  the 
electron  scattering  Is  Isotropic  and  varies  from  it  slightly  if  it  Is  aniso- 
tropic. The  factor  1 - kTQ/E  is  a scaling  factor  to  ensure  that  no  energy  Is 
lost  for  electrons  with  energy  kTQ.  For  a steady  state  then, 

A(£,E,T0)  - B(E,Tq)  (20) 


so  that 


E2 


■ B(E’V 


m 

a 2 
He 


hw 


(21) 


This  relationship,  with  appropriate  values  of  t(E)  and  hw  is  depicted  graphi- 
cally in  figure  15.  Note  In  figure  15,  for  low  values  of  electric  field 
equation  (21)  Is  satisfied  for  two  values  of  the  electron  energy.  For  field 
Ej  consider  an  electron  with  energy  E such  that  £ < £r  Under  the  action  of 
the  field  the  energy  would  increase  to  E Similarly,  for  £x  < £ < Es  the 
losses  predominate  until  the  electron  energy  is  reduced  to  Er  For  fields 
greater  than  EH,  A Is  always  greater  than  B and  the  energy  of  the  electron  will 
Increase  continuously  under  the  applied  field.  This  view,  which  treats  elec- 
trons whose  energy  does  not  lie  to  the  right  of  the  Intersection  of  curves  A 
and  B,  Is  referred  to  as  the  von  Hlppel  low  energy  criteria  and  the  field 
Is  the  breakdown  field.  However,  Frdllch  (ref.  34)  raised  objection  to  this 
treatment  because  It  falls  to  properly  account  for  the  possibility  of  Initially 
present  high  energy  electrons. 

Frdllch  argues  that  when  E > E2  for  E • E2,  then  A Is  greater  than  B and 
this  high  energy  electron  will  continuously  gain  energy  from  the  field.  This 
being  the  case,  £H  has  no  significance  as  a breakdown  field.  To  prevent  such  a 
runaway  condition,  which  due  to  the  shape  of  electron  distributions  can  occur 
at  any  field  no  matter  how  small,  additional  loss  mechanisms  must  be  Introduced. 
The  proponents  of  collective  breakdown  Introduce  free-free  electron  collisions, 
while  proponents  of  avalanche  breakdown  Introduce  free-bound  electron  colli- 
sions. 

Stratton  (ref.  35)  discussed  collective  breakdown  by  observing  that  at 
sufficiently  high  free  electron  densities,  free-free  collisions  will  have  a 
moderating  effect  on  electrons  which  have  energies  which  exceed  the  Frtillch 
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objection.  He  states  that  if  the  free- free  collisions  exceed  the  occurrence  of 
electron-phonon  collisions,  the  electron  distribution  in  the  steady  state  will 
be  Maxwellian  with  a temperature  T greater  than  the  lattice  temperature  TQ. 
Equation  (20)  Is  thus  replaced  by 


*(E.T.T0)  - ff(T,T0) 


where  the  bar  denotes  average  values.  £ and  7 behave  qualitatively  the  same  as 
A and  8 so  that  there  will  be  a critical  field  Ec  which  corresponds  to  the 


electron  temperature  T, 


Is  similar  to.  though  a bit  less  than  von  Hlppel's 


Eh.  The  theory  Is  only  applicable  for  electron  densities  sufficiently  great  so 
that  the  rate  of  energy  loss  by  the  electron  to  the  lattice  Is  less  than  the 
rate  of  energy  loss  by  free- free  Interactions  given  by  Pines  (ref.  36)  as 

B a AffNo** 

e p (23) 

where  q*  Is  an  effective  charge  such  that  (q*/qe)2  * lVe" where  e Is  a dielec- 
tric constant  which  depends  on  the  velocity  of  the  electron.  Thus  the  electron 
density  for  collective  breakdown  must  exceed 


Nc  5 B(E*T0>  = 


Stratton  further  shows  that  the  critical  electron  density  Is  of  the  order  of 
10li  cnf  3.+ 

It  has  been  well  established  In  gas  breakdown  (especially  In  the  photographs 
of  Raether)  (ref.  37)  that  Ionization  proceeds  by  accelerating  electrons  to 
Ionizing  potentials  so  that  collisions  with  neutral  atoms  frees  an  electron, 
which  provides  additional  low  energy  electrons  to  start  the  process  over.  Thus 
Frdllch  assumed  that  for  densities  below  Nc,  electron-neutral  atom  Ionizing 
collisions  would  be  Important  In  solids.  Electrons  of  energy  Ea  greater  than 
the  Ionizing  energy.  I.  (see  figure  15)  can  Ionize  lattice  atoms  and  produce 


This  particular  number  appears  In  many  of  the  theoretical  and  experimental 
studies  of  dc  and  laser- Induced  breakdown. 
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two  low  energy  electrons  per  Ionizing  collisions.  In  steady  state  this  process 
Is  balanced  by  the  reverse  process  of  recombination  In  which  two  free  electrons 
produce  a bound  electron  and  a high  energy  electron.  For  electrons  whose 
energy  satisfies  E2  < E < I they  can  only  increase  In  energy  until  they  reach 
or  surpass  I.  Thus  a steady  state  is  realized  only  when  E2  » I and  equation 
(21)  then  becomes 


E 


2 

I 


«». y 

“Tnr 


q7TU)Tt(ni2Nu+l) 


(25) 


This  Is  referred  to  as  the  Frill Ich  high  energy  criterion.  Stratton  states  that 
once  electrons  achieve  energies  In  excess  of  Ea,  they  are  rapidly  accelerated 
(A  > B)  up  to  energies  of  the  order  I.  Ej  is  then  a lower  limit  for  fields 
which  give  appreciable  Ionization  while  the  actual  breakdown  field  (>Ej)  Is 
determined  by  the  amount  of  Ionization  which  Is  required  to  give  avalanche 
breakdown. 


In  an  electron  avalanche  In  the  simplest  medium  (l.e. » gas)  for  a given 
field  E,  one  electron  will  produce  a electrons  by  collislonal  Ionization  per 
centimeter  of  travel.  The  Townsend  Ionization  coefficient,  o,  is  related  to 
the  reduced  field  In  the  gap.  E/p,  where  p Is  the  pressure  as 


a/p  - (E/p)m 


(26) 


where  m Is  material  dependent.  It  has  a value  of  9.2  for  air  as  reported  by 
Meek  and  Craggs  (ref.  38).  The  growth  of  Ionization  Is  governed  by 


Jx  mHa 

which  for  an  Initial  concentration  of  free  electrons  N0  gives  the  number  N 
after  the  avalanche  has  traveled  a distance  d as 


(27) 


N - N0e 


ad 


(28) 


Breakdown  occurs  when  N reaches  a critical  value,  Nc,  and  from  the  form  of 
equation  (26)  It  Is  easy  to  see  why  a small  increase  In  the  electric  field 
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near  the  breakdown  field  leads  to  a rapid  onset  of  breakdown.  Much  the  same 
view  cf  the  avalanche  process  Is  held  by  those  who  support  the  avalanche  theory 
of  dc  breakdown  In  solid  dielectrics.  The  flrct  problem  encountered  Is  that 
there  does  not  exist  a Townsend  Ionization  coefficient  for  solids  which  exhibits 
the  simple  form  of  equation  (20)  but  rather  one  with  a very  complicated  depen- 
dence on  electric  field  which  cannot  be  reduced  to  quadratures.  For  solids 
Seitz  (ref.  33)  has  given  an  estimate  of  the  critical  electron  density  necessary 
to  produce  a discharge.  He  argues  that  In  a field  of  1 MV/cm  an  electron 
drifting  with  the  field  1 centimeter  while  executing  random  thermal  motion  about 
the  direction  of  drift  sweeps  out  a volume  of  s 3 x 10“*  cm1  and  that  It  seems 
safe  to  assume  that  breakdown  will  occur  If  each  primary  electron  produces  one 
electron  for  each  10s  atoms  In  the  drift  tube.  Thus  when  an  electron  density 
of  * 10* * cm*1  Is  reached,  breakdown  will  occur.  The  total  number  of  electrons 
In  this  example  Is  s 1012  and  since  each  Ionizing  collision  produces  two  elec- 
trons, It  takes  1 generations  to  produce  21  ■ 1012  electrons.  Since  1 here  must 
be  40,  this  Is  the  basis  of  the  Seitz  40  generation  model. 

The  mean  energy  of  an  electron  (figure  15)  must  be  raised  to  I before  an 
Ionizing  collision  can  occur.  Under  the  action  of  the  field  If  the  electron 
experiences  no  collisions.  It  would  be  accelerated  to  an  energy  I In  time 

m(Vr  - vt) 

*1  5E (29) 


or  after  moving  a distance 


ir 


(30) 


In  the  direction  of  the  field.  Typical  values  of  6 eV  for  Ionization  energy 
and  1.0  MV/cm  In  equation  (29)  yield  an  ionization  time  of  about  8 x 10“lH  sec 
or  about  23  full  cycles  of  1.06-um  radiation. 

If  °(tj)  the  probability  that  the  electron  makes  no  collision  for  a time 
tj  „r.d  Tt(£,;  is  the  mean  time  between  collisions,  then  the  mean  rate  of 
Ionization  per  unit  time  Is 

HU) 


45 


the  Ionization  rate  par  unit  path  length  It  thus 


w p(*i> 

° * wjz  ’ Tt(Txm;) i 


(31) 


Breakdown  then  occurs  when  aL  ■ 1 where  L 1$  the  thickness  of  the  sample.  Upon 
using  equation  (31), 


P(tj)L 

1 * Tt(rju(^E 


(32) 


A formula  for  P(t^)  can  be  obtained  by  calculating  the  probability  that  there 
Is  no  collision  between  t and  t + dt. 


P(t  + dt)  - P(t)(l  - dt/t) 


or 

P(tj)  « exp  (-H/E) 


(33) 


(34) 


(35) 


where  H Is  a constant  field.  The  breakdown  field  can  then  be  determined  as 

E»  - nnT^F^rE.Tr  (36) 

Here  one  observes  that  the  breakdown  field  depends  on  sample  thickness  while  In 
collective  breakdown  It  does  net. 

Thermal  breakdown  theories  share  the  final  stage  of  each  of  the  previously 
described  processes.  Obviously,  the  existence  of  an  electron  density  of  101# 
car*  Is  not  reason  In  Itself  for  damage  In  a solid.  Good  conductors  routinely 
have  elactron  dent  1 tits  at  room  t«a*erature  Several  order!  of  magnitude  above 
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1019  air3  with  no  detrimental  effects  (ref.  39).  Rather  It  Is  the  Increased 
conductivity  of  the  dielectric  which  couples  with  the  breakdown  field  (~  1 
MV/cm)  through  joule  heating  to  locally  melt  the  dielectric.  Thermal  breakdown 
theories  assume  that  the  small  conductivity  of  a dielectric  allows  some  joule 
heating  at  the  large  fields  of  Interest  In  breakdown  studies.  The  joule  heating 
raises  the  lattice  temperature  which  raises  the  temperature  dependent  conduc- 
tivity (probably  through  Increases  in  the  number  of  electrons  In  the  conduction 
band  as  well  as  Increases  In  mobility)  and  this  spiraling  process  culminates 
In  localized  melting  and  dielectric  failure. 

Any  calculation  of  the  thermal  breakdown  field  starts  with  the  basic  heat 
flow  equation 


Cy  ■ oE2  - div(<  grad  TQ)  (37) 

but  Its  solution  Is  very  difficult  since  TQ  Is  a function  of  position  and  the 
electrical  conductivity  a .Is  a complicated  function  of  TQ.  k Is  the  thermal 
conductivity.  Usual  attempts  at  solution  are  the  steady  state  In  which  dTQ/dt 
* 0 or  the  short  time  solution  In  which  thermal  conduction  Is  ignored. 

The  conclusions  of  thermal  breakdown  theories  are  generally  valid  at  high 
temperatures  where  the  ratio  of  electrical  to  thermal  conductivity  Is  large. 

This  is  also  true  for  large  electron  densities  such  as  might  be  generated  by  any 
of  the  other  competing  mechanisms. 

Before  considering  theories  of  laser- induced  breakdown  a few  thoughts  drawn 
from  Stratton  (ref.  35)  seem  appropriate.  He  states  that  the  collective  break- 
down theory  of  Frfillch  and  Paranjape  (ref.  35)  provides  a straightforward  method 
for  calculating  a well  defined  critical  dielectric  breakdown  field.  Their 
theory  Is  free  from  objections  If  the  initial  assumptions  of  a Maxwellian  elec- 
tron energy  distribution  Is  valid.  Calculations  based  on  collective  breakdown 
strongly  suggest  the  theory  Is  valid  for  alkali  halides. 

Stratton  further  states  that  the  theory  of  avalanche  breakdown  Is  based  on 
so  many  assumptions,  some  of  which  are  difficult  to  justify,  that  no  great  con- 
fidence can  be  attached  to  the  results.  He  finally  states  that  It  is  possible 
that  avalanche  multiplication  In  dielectrics,  with  a low  Initial  electron  den- 
sity, will  produce  electron  densities  In  excess  of  the  critical  density  Nc  for 
collective  breakdown  before  avalanche  breakdown  takes  place. 
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2.  LASER- INDUCED  BREAKDOWN 

Theories  of  User-Induced  breakdown  of  solid  dielectrics  have  attempted  to 
explain:  (1)  the  mechanism  for  Initiating  primary  electrons,  (2)  the  mechanisms 

responsible  for  creating  electron  densities  of  sufficient  magnitude  to  cause 
damage,  (3)  the  close  correlation  between  dc  breakdown  results  and  laser- 
induced  results,  (4)  the  pulse  width  dependence  of  laser-induced  damage,  (5)  the 
statistical  fluctuations  In  the  breakdown  field  values,  and  (6)  the  frequency 
dependence  of  laser-induced  damage.  Proposed  mechanisms  Include  Inverse 
bremsstrahlung,  collective  Ion  oscillation,  and  various  multiphoton  schemes. 
Basically,  the  current  status  of  the  theory  of  laser-induced  breakdown  Involves 
those  who  confirm  electron  avalanche  as  the  breakdown  mechanism  and  those  who 
claim  that  electron  avalanche  Is  too  Improbable  to  be  operative. 

Bloembergen  (ref.  40)  gives  the  rate  of  increase  in  the  electron  density  N 
during  low-frequency  breakdown  (1 .e. , below  a frequency  at  which  multiphoton 
processes  are  probable)  as 


9N 

at 


a(E)N  + 


tunnel 


(38) 


The  first  term  on  the  right  side  of  the  equation  represents  the  familiar 
avalanche  multiplication  where  a(E)  Is  the  probability  per  unit  time  that  an 
electron  In  field  E will  have  an  Ionizing  collision.  For  the  high  frequency 
limit  the  tunneling  term  becomes  the  same  as  the  multiphoton  Ionization  process 
(ref.  40).  It  Is  this  term  which  Bloembergen  claims,  In  the  absence  of  Initial 
electrons.  Initiates  the  breakdown  process.  Neglecting  the  final  term  In  the 
short  pulse  lengths  of  Interest  and  assuming  an  Initial  electron  density  NQ  of 
10*  - 10l#  cm’1,  one  can  neglect  the  second  term  and  the  result  Is  the  familiar 


N • N. 
o 


W 


(39) 


Bloembergen  then  states  that  at  electron  densities  greater  than  101'  cm*1  the 
energy  deposition  rate  by  absorption  of  the  laser  beam  becomes  so  high  that  the 
temperature  rise  of  the  lattice  Is  sufficient  to  cause  damage  even  In  picosecond 
pulses. 
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Fradin  et  al.,  (ref.  17)  further  refined  B1 oembergen ' s argument  by  stating 
that  the  process  of  electron  avalanche  involves  energy  exchange  between  the 
field  and  electrons  which  is  governed  by  the  formula  for  ac  conductivity 


d£  . "VT 

^ m(1+w*Ta) 


and,  since  the  details  of  energy  input  determine  the  energy  distribution  and 
thence  N,  that  threshold  field  scales  with  frequency  as 


El  - 1 ♦ uiV  (41) 

where  t Is  the  electron-phonon  collision  relaxation  time.  Equation  (41)  is  used 
to  explain  the  observed  factor  of  1.5  between  dc  and  optical  frequency  breakdown 
fields. 

Bass  and  Barrett  (ref.  41)  propose  a probaballstlc  view  of  laser.;- induced 
damage  based  on  (cf.  equation  (35)) 


Pt  « exp  (-K/E)  (42) 

as  the  probability  that  a single  pulse  will  produce  damage  at  the  given  field  E. 
The  avalanche  theory  Is  modified  to  account  for  the  possibility  that  an  electron 
will  suffer  only  lucky  collisions;  that  is,  collisions  which  reverse  the  momen- 
tum of  the  electron  just  as  the  electric  field  reverses.  Such  collisions  are 
necessary  to  Increase  the  energy  of  a free  electron  In  an  alternating  field  since 
the  momentum  of  the  electron  In  an  alternating  field  and  the  phase  of  the  elec- 
tric field  are  always  opposed. 

The  pulse  length  dependence  was  treated  by  Fradin  (ref.  16)  using  equation 
(39).  Assuming  an  Initial  electron  density  of  10*.requ1red  a multiplication 
Mc  by  avalanche  of  10xo.  Then  the  Ionization  constant  Is  taken  as 

o(E)  - *nMc  - 18/tp  (43) 

where  tp  is  the  pulse  duration. 
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Without  belaboring  the  subject,  this  review  will  be  concluded  by  citing 
the  highly  readable  paper  by  Sparks  (ref.  42}  on  his  objections  to  the  current 
avalanche  theories  of  laser-induced  breakdown.  Sparks  shows  that  the  avalanche 
multiplication  constant  Is  too  small  by  tens  of  orders  of  magnitudes  In  some 
cases  to  explain  damage  results  even  when  Initial  electron  densities  as  high  as 
10* 0 cm-3  are  assumed.  The  proposed  frequency  dependence  of  the  breakdown  field 
E£  ~ 1+o)2t2  disagrees  with  experimental  results.  The  explanation  In  terms  of 
an  anomalously  small  electron  relaxation  time  was  shown  to  be  Inconsistent  with 
the  value  of  t required  to  explain  the  magnitude  of  Eg  with  calculated  values 
of  t and  with  the  difference  between  the  dc  and  1.06-ym  breakdown  results. 

Sparks  further  shows  that  the  temperature  dependence  of  Eg  Is  Incorrect  and 
that  values  of  NQ  In  the  range  10*-1010  cm"3  would  place  the  alkali  halides  In 
the  semiconductor  class  of  materials.  Further,  densities  of  this  magnitude  are 
Inconsistent  with  photoconductivity  measurements  and,  in  fact,  limits  on  NQ  set 
by  conductivity  arguments  estimate  the  probability  of  finding  even  one  electron 
In  the  focal  volume  of  typical  laser  experiments  as  being  less  than  10”*.  The 
lucky  electron  theory  was  considered  and  even  with  overgenerous  estimates  of 
the  probability  of  backscatterlng  Into  phase  with  the  field  yielded  a probabil- 
ity of  the  success  of  an  electron  reaching  an  Ionizing  energy  of  ~10"52.  Sparks 
then  proposes  a preliminary  new  theory  Involving  Holstein  processes  (photon- 
electron-phonon)  and  avalanche  mechanisms. 

With  these  remarks  as  background,  a new  approach  to  laser- Induced  damage 
utilizing  direct-field  Ionization  will  be  proposed. 

3.  AN  ALTERNATIVE  TO  AVALANCHE  BREAKDOWN 

Because  of  the  objections  to  an  avalanche  theory  of  solid  dielectric  break- 
down voiced  by  Stratton  and  the  equally  vigorous  objections  to  an  avalanche 
theory  of  laser- Induced  damage  raised  by  Sparks  It  does  not  appear  fruitful  to 
give  here  yet  another  attempt  at  a modification  of  basic  avalanche  theory. 

Rather  a simpler  theory  of  direct-field  Ionization  will  be  proposed  which  has 
the  advantage  over  the  more  complicated  theory  In  that  it  is  predictive.  The 
model  and  basic  assumptions  will  provide  a basis  for  the  mathematical  treatment 
to  follow. 

a.  Assumptions 

One  of  the  few  threads  cownon  to  nearly  all  theories  of  dc  breakdown 
and  laser-induced  damage  to  solid  dielectrics  Is  that  an  electron  density  of  the 
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order  of  1018  cm-*  Is  required  before  energy  conversion  to  the  lattice  through 
absorption  by  the  electrons  Is  sufficient  to  damage  the  lattice.  David  et  al., 
(ref.  43)  measured  electron  densities  of  ~!0ia  cm*5  required  to  effective  Ijrtfc- 
gln  to  absorb  laser  energy.  Guenther  and  Pendleton  (ref.  44)  reported  densities 
of  1018  cm*3  at  visible  plasma  thresholds  In  deuterium  gas.  Recently  Alyasslnl 
and  Parks  (ref.  45)  calculated  the  electron  density  at  surface  damage  threshold 
using  the  observed  change  In  refractive  Index.  The  values  they  report  are 
~1018  cm*3.  Thus  a basic  assumption  Is  that  If  electron  densities  of  the 
order  of  1018  cm’3  can  be  provided  by  the  laser  pulse,  damage  will  follow  by 
the  subsequent  absorption  of  the  remainder  of  the  pulse. 

The  other  basic  assumption  around  which  the  theory  Is  devised  is  that  a 
semlclasslcal  description  of  the  field-atom  Interactions  can  provide  the  req- 
uisite electron  density. 

b.  The  Model 

The  author  proposes  here  a model  which  explains  the  salient  features  of 
pulsed,  laser-induced  damage  to  transparent^ dielectrics.  The  laser  beam  liber- 
ates electrons  from  the  atoms  In  th^  region  of  Interaction  between  the  laser 
beam  and  the  material.  As  has  been  stated  previously,  It  Is  not  the  liberation 
of  electrons  which  damages  the  dielectric  but  rather  the  deposition  of  the 
energy  of  the  laser  beam  In  the  electronic  plasma  which  heats  the  lattice  and 
causes  damage.  The  behavior  observed  by  Milam,  et  al.,  (ref.  45)  supports  this 
conjecture.  Milam  determined  the  starting  times  of  damage  events  by  observing 
the  attenuation  of  the  laser  pulse  which  was  transmitted  by  the  dielectric.  The 
majority  of  starting  times  occurred  at  or  before  the  peak  laser  Intensity.  A 
few  damage  events  occurred  after  the  peak  of  the  laser  pulse,  but  In  no  case  did 
damage  occur  In  the  final  third  of  the  laser  pulse  duration.  One  explanation 
of  these  observations  Is  that  once  a critical  electron  number  density  is  reached, 
sufficient  to  Initiate  absorption  of  the  laser  beam,  damage  will  occur  only  If 
a significant  portion  of  the  energy  In  the  laser  beam  is  available  for  absorption 
and  subsequent  heating  of  the  lattice  to  localized  destruction.  A heuristic 
argument  can  be  given  here  for  the  observed  Increase  In  threshold  damage  fields 
with  decreasing  laser  pulse  length  (ref.  47).  If  one  assumes  that  a constant 
heat  Input  Is  necessary  to  Induce  damage  In  a given  material,  then  the  heating 
rate  must  vary  Inversely  with  the  duration  of  the  pulse.  Since  the  heating 
rate  Is  proportional  to  the  Input  power  density,  the  power  density  and  thus  the 
threshold  field  will  also  exhibit  an  Inverse  pulse  length  dependence.  It  Is  also 
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possible  to  suggest  that  there  is  a maximum  pulse  length  for  which  the  Inverse 
pulse  length  dependence  is  valid.  For  pulse  lengths  greater  than  electron 
recombination  times,  an  increase  in  the  energy  of  the  laser  pulse  is  necessary 
to  overcome  the  'losses  In  the  electron  density  due  to  electron  recombinations. 
Thus  one  would  expect,  if  not  a leveling  off,  at  least  a decrease  In  the  inverse 
pulse  length  dependence. 

Finally,  although  the  tacit  assumption  has  been  made  that  linear 
absorption  can  be  neglected,  the  proposed  model  can  be  used  to  indicate  the 
expected  behavior  of  threshold  fields  on  linear  absorption.  Quite  simply,  the 
linear  absorption  of  the  laser  beam  by  the  dielectric  enters  directly  into  the 
heating  rate  of  the  lattice.  Thus  one  would  expect  that  an  Increase  in  linear 
absorption  woul.d  result  in  a decrease  in  the  threshold  field  for  laser-induced 
damage. 

c.  Alternative  Theory 

Consider  a collection  of  N atoms  per  cubic  centimeter.  Let  an  electric 
field  E - E0eJa,t  be  incident  upon  the  material.  Further  consider  the  inter- 
action of  this  field  with  an  average  optical  electron  on  a typical  atom  in  the 
collection.  The  motion,  in  the  classical  charged  particle  on  a spring  model, 
of  the  electron  is  given  by  solution  of  the  phenomenological  equation  of 
motion 


m{x*+  yk  + u*x)  a F * (44) 

where  y is  an  effective  damping  or  friction  coefficient,  u)Q  is  the  natural 
oscillation  frequency,  x is  the  displacement,  the  dot  denotes  differentiation 
with  respect  to  time,  F is  the  applied  force,  and  tha  subscript  A denotes  local 
field  at  the  atom.  Since  E varies  sinusoidally  the  displacement  will  oscillate 
at  the  driven  frequency  and  substituting  x a xQeJwt  into  equation  (35)  yields 

qe/m 

X * -w^  + jytL)  + ^A  (45) 

The  induced  dipole  moment  of  an  atom  p is  qax  thus 

* ^ . 

p * Jvui  + 1*7^  l (46) 
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For  linear,  isotropic  media  the  polarization  is  proportional  to  the  applied 
field  with  constant  of  proportionality  ct(w),  the  atomic  polarizability.  Thus 
we  write 


eoo(“)  *4 


(47) 


Then  using  equation  (46)  we  cbtain 


a(w) 


qe2/m  eO 

-ui2  + jycj  + oj2 


(48) 


This  simple  expression  is  quite  close  to  the  rigorous  expression  obtained  from 
quantum  mechanical  considerations  and,  in  fact,  can  be  transformed  by  only  a 
few  modifications  into  the  more  correct  form.  The  differences  In  the  classical 
approach  and  the  quantum  mechanical  are  that  each  atom  has  not  one  natural  fre- 
quency of  oscillation  but  several  natural  frequencies  each  of  which  responds 
to  the  field  with  a different  damping  coefficient  yk  and  a different  oscillator 
strength  fk,  Equation  (48)  is  transformed  by  these  notices  into 


a(ui) 


3L  It 

V k + * “k 


The  volume  polarizability  is 


(49) 


P - M«)*0Et  (50) 

and  was  used  in  the  formulation  of  Lorentz  as  reported  by  Jackson  (ref.  48) 
to  qbtain  the  effective  local  field,  E^,  of  an  atom  imbedded  in  an  Isotropic 
medium  where  a macroscopic  field  E exists.  For  fields  whose  wavelength  is  much 
longer  than  the  spacing  between  atoms  the  local  field  of  the  atom  in  an  iso- 
tropic medium  is 


E 


0 

•-* 


■ E + 


P 


(51) 
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Combining  equations  (50)  and  (51)  results  in  the  Clausius-Mosotti  or  Lorentz- 
Lorenz  formula 


n2-l 

n*+F 


Na(u>) 


(52) 


where  n is  the  macros copically  averaged  quantity,  the  refractive  index. 
Now  equation  (45)  can  be  rewritten  as 


£oc<(ai)  h 


(53) 


or  by  using  equation  (52)  this  becomes 


nM  3(£o> 


n*+2  Nqe  "l 


(54) 


as  the  relationship  between  the  disp^cement  of  an  electron  and  the  local  field 

Er 

Rewriting  this  expression  in  terms  of  the  macroscopic  field  by  the 
Lorentz  local  field  correction  yields 


E * -JL  !!§.  x 

"a-l  «o 


(55) 


To  obtain  the  threshold  for  laser-induced  damage,  one  must  account  for  the  fact 
that  a critical  free  electron  density  of  Nc  cnr3  must  be  obtained  by  the  field 
acting  on  the  atoms  and  one  must  decide  a reasonable  critical  displacement, 
xcr  bey°nd  wh,*cb  *be  electron  is  "freed"  from  the  atom. 

If  a given  total  energy  is  required  to  iorlze  every  atom  In  the  focal 
region  then  only  a portion  Nc/N(cnr3)  will  be  required  to  ionize  Nc  cm”3. 

Since  the  optical  electric  field  scales  as  the  square  root  of  the  total  energy, 
equation  (55)  becomes 


*th 


N qe 

ri*rr  IT  xcr  1 r* 


(56) 
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Physical  arguments  can  now  be  used  to  obtain  an  order  of  magnitude 
estimate  for  the  critical  displacement  xcr.  A lower  limit  can  be  obtained  from 
the  arguments  of  Wunsch  (ref.  49).  He  shews  that  the  electric  dipole  moment  of 
an  atom  under  the  action  of  an  applied  field  does  not  saturate  until  the  elec- 

O 

tron  Is  displaced  0.75  - 1.0  A for  most  optical  materials.  Physically  this 
Implies  that  the  restoring  force  assumed  in  deriving  equation  (55)  through  the 
electric  dipole  Interaction  lis  valid  until  saturation  occurs  and  thus  that  the 
electron  Is  still  Intimately  attached  to  the  atom.  An  upper  limit  for  xcr  can 
be  obtained  by  observing  that  on  the  average  an  atom  In  a solid  Is  separated 
from  Its  nearest  neighbor  by  a distance 

as  - N-1/3  (57) 

O 

For  ordinary  materials  this  distance  corresponds  to  2 to  3 A.  It  Is  not  there- 
fore unreasonable  to  assume  that  when  an  electron  Is  half  way  between  two  atoms 
it  Is  essentially  liberated.  Thus  the  critical  distance  becomes 

x » i 

xcr  T (58) 


which  Is  consistent  with  the  lower  limit  described  above. 

The  final  expression  for  the  threshold  electric  field  In  terms  of 
material  properties  then  becomes 


N % N2ZiN-!/2  .VNc  rrr3  V/m 
n2-l  eQ  2 


(59) 


or  combining  terms 


Eth  * n^T  ^ WV"*-3  v/m  (60) 

It  will  prove  Instructive  later  to  write  equation  (60)  In  terms  of  the  average 
separation  as 


(61) 
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To  obtain  an  order  of  magnitude  value  for  the  predicted  threshold  field  the 
quantity  1018  cm"3  (1024  m"3)  Is  substituted  for  the  critical  electron  density 
N In  equation  (61)  to  ootain 

V 

E . 3,-°5.x  10 3 MV/Cffl  (62) 

th  ("*-’) 

As  will  be  shown  below,  equation  (62)  adequately  models  the  behavior 
of  a variety  of  materials  subjected  to  damaging  laser  irradiation.  Equation 
(62),  however,  was  developed  to  explain  laser-induced  damage  in  the  bulk  of 
materials.  Because  the  experimental  portion  of  this  study  dealt  with  surfaces 
of  materials,  additional  considerations  must  be  taken  intj  account  to  make 
equation  (62)  directly  applicable  to  this  report. 

Surfaces  introduce  additional  features  such  as  contamination,  surface 
finish,  and  roughness.  A basic  assumption  In  applying  equation  (62)  to  surfaces 
Is  that  the  process  of  laser-induced  damage  is  the  same  for  surface  and  bulk 
regions  of  dielectrics.  For  a perfect  surface  the  average  atomic  spacing  and 
the  surface  roughness  should  correspond  directly.  In  such  an  event,  one  would 
expect  that  the  surface  and  bulk  threshold  fields  would  be  the  same.  Thus, 
equation  (62)  Is  directly  applicable  for  perfect  surfaces.  Equation  (62)  Is 
still  of  use— lacking  perfection  but  possessing  similarity.  If,  Indeed,  the 
process  of  laser-induced  damage  is  the  same  for  surfaces  and  bulk,  then  a set 
of  different  materials  with  Identical  physically  realizable  surface  finishes 
(in  terms  of  roughness  and  freedom  from  contamination)  will  still  scale  in 
damage  threshold  with  equation  (62).  One  would  expect  Imperfect  surfaces  to 
exhibit  damage  thresholds  less  than  those  of  the  corresponding  bulk  materials, 
but  the  relative  values  between  materials  should  still  hold.  The  problem 
encountered  Is  one  of  a lack  of  similarity.  That  Is,  the  surface  of  the 
several  materials  tested  in  this  study  were  not  Identical. 

In  a study  of  laser-induced  damage  as  a function  of  surface  roughness 
House  (ref.  50)  showed  that  for  a given  material  the  threshold  electric  field 
for  laser-induced  damage  scaled  with  surface  roughness,  o,  as 

V®  Eth  * constant  (63) 
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Then  for  two  samples  with  surfaces  roughness  and  o2,  the  threshold  fields 
E!  and  E2  are  related  by 


(64) 


House  found  that  this  relationship  held  for  the  ratio  of  the  surface  threshold 
of  a material  with  roughness  a to  the  bulk  threshold  of  that  material  with 
average  atomic  spacing  o$  as 


^surface 


Ebulk 


(65) 


Then  to  compare  materials  with  different  surface  roughnesses,  one  chooses  a 
standard  roughness  value,  and  multiplies  the  measured  threshold  field 
by  the  square  root  of  the  roughness  ratio 


Vf^  E * ESTAND  (66) 

where  E,  a are  the  measured  threshold  field  and  surface  roughness,  and  657/^0 
Is  the  threshold  field  for  a sample  with  the  standard  roughness.  Thus  the 
material -to-materlal  variation  of  the  threshold  electric  field  can  be  Isolated 
by  comparing  the  various  materials  at  an  equivalent  standard  roughness. 

As  a final  topic  consider  the  pulse  length,  spot  size,  and  laser  fre- 
quency dependence  of  laser-induced  damage.  These  are  contained  In  the  details 
of  the  plasma  heating  by  the  laser  pulse.  To  calculate  the  heating  rate  of  a 
plasma,  suppose  that  laser  power  density  S Is  Incident  upon  a plasma  of  thick- 
ness dx  with  absorption  coefficient  K.  Then  by  Beer's  Law  the  power  density 
after  traversing  the  plasma  is  S e-Kdx.  The  absorbed  energy  then  Is 

dE  ■ ASAdt  ■ (S  - Se-Kdx)Adt  (67) 

where  A Is  the  area  of  the  laser  beam,  AS  Is  the  change  In  S In  traversing  a 
distance  dx  of  the  plasma,  and  dt  Is  a differential  elements?  time.  Expanding 
the  exponential  for  small  Kdx  yields 
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dE  • 5[1  - (1  ♦ Kdx  + ...))  Adt 


dE  = KS  dx  Adt 


Then  the  heating  rate  Is  given  by 


, dE 
* 3t 


Adz  K S 


A given  material  will  fall  when  the  total  heat  Input  to  the  lattice  ;>  suffi- 
cient to  cause  localized  melting.  For  a given  material  then  the  total  heating 
is  a constant  and  equation  (70)  yields 


constant 


tp  Adx  K S 


If  the  thickness  of  the  plasma,  dx.  Is  constant  at  threshold  for  a specific 
material,  then  equation  (71)  can  be  used  to  Investigate  the  pulse  length  and 
spot  size  dependence  of  laser-induced. damage.  The  absorption  In  electronic 
plasmas  Is  governed  by  the  free-free  or  Inverse  Bremsstrahlung  process  which 
has  an  absorption  coefficient  given  by  Johnston  and  Dawson  (ref.  51)  of 


3.08  x 1C*7  Z N*  tnA(w) 

K . i- 

(KT  )•/» 


« 

«p 


where  kT#  is  In  eV,  A (w)  • Vy/aip  p(rt(|.  Here,  Zqt  Is  the  Ionic  charge,  T#  Is  the 
electron  temperature,  N#  1$  the  electron  number  density*  VT  Is  the  thermal 
velocity  of  the  electrons,  pm1n  Is  the  minimum  Impact  parameter  for  electron-ion 
collisions  [pm1n  £ maximum  of  ZqJ/kT  or  (h/mkT)1/1). 

To  quantify  equation  (72)  It  was  shown  by  David  et  al.,  (ref.  43)  that 
both  Nq  and*T  obeyed  a power-law  dependence  with  laser  flux.  Faugulgnon  and 
Floux  (ref.  52)  calculated  and  experimentally  verified  that  t#  - $4/#.  Boland 
et  al.,  (raf.  53)  showed  that  the  expanding  plasma  obeyed  ah  adiabatic  expansion 
law  and  verified  experimentally  that 
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Te  ■ const  Nft2/3 


(73) 


Then 


TVT  - S (74) 

e ' 

For  a given  materiel  at  a fixed  laser  frequency  and  constant  spot  size,  equa- 
tions (71)  and  (74)  can  be  used  to  write 


v 


2 


v-;nst 


(75) 


or 


E • const  tp-‘/‘  (76) 

Thus,  the  threshold  electric  field  varies  Inversely  with  the  fourth  root  of  the 
pulse  length. 

For  a given  material  at  a fixed  laser  frequency  and  constant  pulse 
duration  equations  (71)  and  (74)  can  be  used  to  write 

A S*  ■ const  (77) 


or 


E « const 

JS 


(78) 


where  d Is  the  diameter  of  the  focal  spot  of  the  laser.  The  limits  on  the 
applicability  of  equations  (76)  and  (78)  will  now  be  considered. 

There  are  probably  a lower  and  upper  limit  on  the  pulse  duration  for 
which  equation  (76)  Is  expected  to  hold.  Laser  frequencies  are  of  the  order  of 
10u  - 10ls  sec”1.  Such  atomic  processes  as  tunneling  require  about  10"1*  sec. 
Thu:>,  equation  (76)  will  probably  not  hold  for  pulse  durations  less  than  the 
ordor  of  10* u sec.  This  Is  about  two  orders  of  magnitude  less  than  the  shortest 
pulne  durations  currently  available,  and  thus  the  conjecture  cannot  now  be 
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tested.  Because  equation  (76)  was  derived  by  Ignoring  losses,  it  is  clear  that 
loss  mechanisms  will  invalidate  the  predicted  behavior.  A principle  loss 
mechanism  Is  electron  recombination.  Since  recombination  times  are  of  the 
order  of  microseconds,  equation  (76)  Is  not  expected  to  hold  for  pulse  dura- 
tions much  greater  than  a few  microseconds.  Thus,  the  lower  time  limit  would 
appear  to  be  process  dominated  while  tne  upper  time  limit  is  loss  dominated. 

The  spot  size  dependence  will  also  exhibit  a range  of  applicability. 

For  snail  spot  sizes,  losses  due  to  thermal  diffusion  become  Important.  The 
heat  added  to  the  lattice  diffuses  out  of  the  focal  volume  In  a time  on  the 
order  of  td  ■ r2/D,  where  0 is  the  thermal  dlffusivity  and  r is  the  radius  of 
the  focal  spot.  Losses  by  diffusion  then  become  Important  when  the  pulse 
length,  tp,  is  of  the  order  of  the  diffusion  time.  For  a 40-nsec  pulse  dura- 
tion and  a typical  dlffusivity  of  1 to  10  cm2/sec,  this  implies  that  when  the 
spot  diameter  is  less  than  about  1 micrometer  thermal  diffusion  losses  begin 
to  negate  the  dependence  given  in  equation  (76).  Furthermore,  when  the  spot 
size  becomes  large  enough  a minimum  laser  intensity  Is  reached  such  that 
further  decreases  In  the  intensity  are  not  possible  while  still  maintaining  a 
damaging  level  of  irradiation.  Thus,  there  Is  a minimum  laser  intensity  which 
is  necessary  to  produce  a critical  number  density  of  free  electrons  sufficient 
to  Initiate  damage.  For  an  estimate  of  the  maximum  spot  size  for  which  equa- 
tion (78)  holds,  consider  that  an  intensity  level  typical  of  damage  In  the 
absence  of  absorption  (~109  watts/cm2)  is  the  minimum  intensity  level  suffi- 
cient to  produce  about  1018  electrons/cm3.  A level  of  1017  electrons/cm3  is 
probably  too  low  to  Initiate  intrinsic  damage,  thus  an  Increase  In  the  typical 
100  urn  spot  size  by  a factor  of  10  to  1 mm  will  probably  result  in  little  de- 
crease in  threshold  intensity. 
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SECTION  IV 
RESULTS 

Internal  consistency  Is  lacking  In  the  many  published  reports  on  laser- 
induced  damage.  One  of  the  major  causes  Is  tiie  many  ways  In  which  the  raw  data 
collected  by  the  Investigators  Is  reduced  to  the  final  form.  To  obviate  the 
possibility  that  this  report  will  contain  Irreproducible  results,  the  collection 
and  transformation  of  the  raw  data  Into  final  form  will  now  be  discussed. 

1.  RAW  DATA 

For  this  study  the  variables  of  Interest  Included  laser  energy,  laser  pulse 
duration,  focal  spot  size,  and  whether  damage  occurred  on  a given  shot.  The 
first  and  last  of  these  has  been  discussed  above  under  active  diagnostics  and 
will  not  be  rediscussed  here  except  to  note  that  on  each  Irradiation  either  a 
"damage"  or  a "no  damage"  was  recorded  along  with  the  total  energy  Incident 
upon  the  target.  The  focal  spot  size  was  measured  using  the  technique  described 
In  section  II. 2,  Parameterization.  This  parameter  was  rechecked  periodically 
throughout  the  study  and  was  found  consistent  to  within  ±3.5  percent.  The  tem- 
poral shape  of  the  laser  pulse  was  recorded  on  a fast-sweep  oscilloscope,*  and 
the  pulse  duration  was  taken  as  the  time  separation  between  the  trace  at  half 
the  maximum  deflection. 

2.  DATA  REDUCTION:  THE  ENERGY  THRESHOLD 

A method  for  determining  the  energy  which,  for  a given  sample,  can  be  re- 
ported as  the  energy  threshold  was  developed.  It  Is  a method  which  does  not 
rely  on  subjective  judgement,  and  thus  It  must  therefore  remove  one  of  the 
largest  uncertainties  In  experimental  work.  Before  the  final  method  to  be 
described  here  was  decided  on,  a total  of  1438  shots  were  taken  on  20  samples 
to  assist  In  the  determination  of  an  objective  method  for  obtaining  a damaging 
energy  threshold  value  for  each  sample.  It  should  be  mentioned  here  that  there 
Is  a body  of  evidence  collected  by  Bass  and  Barrett  {ref.  41),  principally,  and 
HI 1am  and  co-workers  (ref.  46)  that  suggest  that  a damage  threshold  does  not 
exist  but  that  the  damage  process  Is  entirely  probabalistic.  Their  view  Is 
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not  unrealistic  and  in  fact  the  damage  thresholds  reported  here  should  corre- 
spond roughly  to  their  P(0.5)  point  or  the  energy  at  which  the  probability  is 
0.5  that  damage  will  occur  on  a single  shot.  Figure  16  illustrates  the  method 
used  in  determining  energy  threshold.  Each  shot  Is  represented  by  an  energy 
and  a +1  or  -1  corresponding  to  damage  or  no  damage,  respectively.  The  highest 
energy  event  which  did  not  cause  damage  Is  then  connected  with  a line  to  the 
lowest  energy  event  which  did  cause  damage  and  the  energy  at  which  the  line 
crosses  the  zero  axis  is  defined  to  be  the  energy  threshold.  This  method  re- 
moves subjective  judgement  from  the  report  of  the  damage  threshold.  By  effec- 
tively averaging  the  highest  no-damage  and  the  lowest  damage,  the  effect  of 
either  taken  singly  is  moderated.  There  are  those  who  use  a similar  method  but 
report  only  the  highest  energy  which  did  not  cause  damage  as  the  threshold,  but 
this  is  much  more  sensitive  to  the  unusual  endpoint  than  Is  the  averaging  method 
espoused  here.  The  number  of  shots  whlcn  one  would  need  to  take  to  use  only  the 
highest  no-damage  would  have  to  be  many  times  the  25  shots  on  each  sample  used 
In  this  study.  In  the  1438  shots  in  20  samples  alluded  to  above,  it  became 
apparent  that  150  shots  per  sample  would  not  enable  one  to  select  a damage 
threshold  with  more  accuracy  than  the  25  shots  used  here  and  was  thus  not 
economically  justified.  However,  very  different  results  were  possible  If  less 
than  20  shots  were  used.  In  particular  15  and  12  shots  could  vary  from  the 
results  of  25  to  150  shots  by  25  to  30  percent. 

The  ratio  R was  computed  as  the  lowest  damage  event  divided  by  the  highest 
no-damage  event.  This  ratio  was  used  as  a measure  of  the  uncertainty  In  the 
result  by  computing 


— y- - ■ % uncertainty  (79) 

In  most  cases  the  uncertainty  In  energy  threshold  ranged  from  ±15  percent  to 
±20  percent,  reflecting  the  fact  that  the  energy  at  which  the  sample  will  damage 
on  the  first  shot  every  time  Is  usually  at  least  twice  the  energy  at  which  the 
sample  will  never  damage.  Thus,  unless  one  can  take  the  several  thousand  shots 
per  sample  needed  to  accurately  map  the  probability  to  damage  versus  energy, 
there  does  not  seem  to  be  a better  method  for  determining  a “damage  threshold." 

3.  DATA  REDUCTION:  THE  ELECTRIC  FIELD  THRESHOLD 

In  this  report  the  threshold  values  which  are  listed  are  the  optical  elec- 
tric field  thresholds.  From  equation  (4)  we  see  that  the  average  electric 
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Figure  16.  Method  Used  to  Determine  Energy  Threshold  for  Damage 
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field  Is  related  to  the  average  power  density  which,  for  a Gaussian  temporal 
pulse,  is  given  by 


tp  • area 


(80) 


where  E is  total  energy,  tp  is  the  full  pulse  width  at  half  maximum  intensity 
and  the  area  of  the  focal  spot  was  taken  as  that  lying  within  a radius  at  which 
the  Intensity  fell  to  half  the  maximum  intensity.  For  this  study  the  half-power 
point  had  a diameter  of  104  ym.  Although  the  optical  electric  field  Is  re- 
ported, It  was  the  energy  threshold  which  was  experimentally  determined.  The 
electric  field  was  then  obtained  using  the  average  pulse  width  (measured  for 
the  particular  sample)  to  calculate  the  average  power  density  from  equation  (80) 
and  the  field  from  equation  (4).  The  pulse  length  variation  was  ±5  percent. 
Shots  which  had  overly  long  pulse  widths  were  discounted  and  not  Included  In  the 
data  reduction.  Less  than  2 percent  of  the  total  shots  we^e  retaken  because  of 
temporal  pulse  width  variation.  Variation  of  the  pulse  duration  usually  indi- 
cated a slight  misalignment.  Realignment  brought  pulse  duration  back  to  normal. 

From  the  continuity  of  the  tangential  electric  field  the  expression 


E - 


2n 


TT 


19.4  Vs 


(81) 


gives  the  field  at  the  surface  of  a dielectric  with  refractive  Index  n.  This 
also  Is  the  field  at  the  surface  of  a nonlossy  dielectric  film  placed  on  the 
surface  In  integral  multiples  of  a half  wavelength  optical  thickness.  Except 
In  the  case  In  which  a film  has  a lower  refractive  Index  than  the  substrate,  it 
also  represents  the  peak  electric  field  In  half-wave  films  such  as  tested  In 
this  study.  It  Is  this  calculated  field  which  is  reported  as  the  threshold 
values  In  this  report. 

4.  DATA  REDUCTION:  ACCURACY  OF  RESULTS 

The  measured  energy  has  an  uncertainty  of  ±7  percent  due  to  the  4 percent 
uncertainty  In  the  pyroelectric  detector  and  the  ±5.5  percent  uncertainty  in  the 
calorimeter  used  to  calibrate  the  pyroelectric  detector.  The  average  uncertainty 
In  threshold  energy  was  ±19  percent.  With  a pulse  width  uncertainty  of  5 percent 
and  an  uncertainty  In  the  focal  spot  diameter  of  3.5  percent  the  electric  field 
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was  known  for  a given  shot  to  ±9.5  percent.  Folding  In  the  19  percent  for  the 
threshold  energy  results  In  a total  uncertainty  in  the  threshold  electric  field 
of  ±13.4  percent. 

5.  TEST  OF  THEORY  VERSUS  PREVIOUS  RESULTS 

Sufficient  Information  is  contained  in  the  papers  by  Smith  (ref.  54), 
DeShazer  et  al.,  (refs.  11,  15)  Austin  et  al.,  (ref.  14),  and  Fradin  et  al., 
(ref.  17)  to  enable  one  to  compare  their  experimental  results  against  theoreti- 
cal predictions  from  equations  (55),  (62),  (76),  and  (78).  For  a given  laser 
with  fixed  frequency  and  pulse  width,  equations  (55)  and  (62)  should  give  the 
relative  damage  thresholds  for  various  materials.  It  is  desirable  to  compare 
the  two  theoret.cal  expressions  In  an  attempt  to  identify  the  major  contribution 
to  the  expected  behavior.  In  so  doing,  the  critical  displacement  x in  equa- 

w r 

tlon  (55)  Is  simply  taken  as  a constant  and  the  experimentally  determined 
threshold  fields  are  then  compared  to  N/n2-1.  For  the  full  expression  given  In 
equation  (62),  the  average  atomic  spacing,  a$.  Is  just  the  cube  root  of  the 
Inverse  of  the  atomic  nu/riber  density. 

The  experimental  study  detailed  In  this  report  employed  a laser  with  a 
40-nsec  pulse  width.  However,  It  Is  Interesting  to  compare  the  pulse  length 
dependence  predicted  by  equation  (76)  with  electric  field  thresholds  experimen- 
tally determined  at  five  different  pulse  lengths.  To  this  end,  it  will  be 
demonstrated  that  the  process  of  laser-induced  damage  at  30  psec  Is  essentially 
the  same  as  the  process  at  40-nsec  pulse  durations.  The  validity  of  equation 
(76)  will  then  be  tested  directly  for  NaCi  at  five  pulse  durations  from  15  psec 
to  10  nsec. 

The  spot  size  dependence  of  laser-induced  damage  will  be  considered  as  the 
final  comparison  of  the  theory  with  previously  published  results. 

In  the  tables  and  graphs  which  follow,  the  field  Is  the  threshold  optical 
field  for  laser- Induced  damage. 

First  consider  the  data  gathered  by  DeSh'aze  * and  co-workers.  They  placed 
various  thin  films  on  BSC-2  glass  substrates.  The  data  have  been  collected  in 
table  1.  These  data  are  plotted  In  figures  17  and  18.  The  least  square  fit 
to  the  data  has  a root  mean  square  variation  of  ±20.8  percent  for  compared 
to  the  full  expression  and  ±29,9  percent  when  compared  to  only  N/n2-l.  Note 
that  the  thicker  films  fall  at  lower  fields,  a result  attributable  to  slight 


65 


AFWL-l'fl-76-61 


absorption  ir.  the  thin  films  according  to  an  industry  source  (ref.  29)  but 
absorption  which  has  not  been  quantified. 

Table  1 


DAMAGE  THRESHOLDS  FOR  FILM 
(adapted  from  reference  11) 


Material 

Thickness 

Index 

Field 
(MV/ cm) 

N 

n^T 

(1028  irf3) 

9.05  x 103 

/cr3(nM) 

(MV/cm) 

Zns 

A/4 

2.35 

0.208 

1.12 

1.22 

T102 

A/4 

2.26 

0.513 

2.0? 

1.43 

ZrO. 

A/4 

1.975 

0.569 

2.78 

2.05 

Si02 

A/4 

1.449 

1.656 

6.02 

5.23 

Ti0. 

A/2 

2. 28 

0.378 

2.07 

1.43 

Zr02 

A/2 

1.975 

0.454 

2.78 

2.05 

SiO, 

A/2 

1.44b 

1.54 

6.02 

5.23 

Ti02 

3A/4 

2.28 

0.369 

2.07 

1.43 

The  data  collected  by  Putman  (ref.  55)  and  presented  by  Austin  et  al., 

(ref.  14)  are  presented  next.  The  material  consisted  of  half  wavelength  (at 

i 

1.06  urr»)  thick  films  comprised  of  vapor-phase  mixtures  of  ThF4  and  ZnS.  The  da 
are  collected  in  tabl*  2 indicating  tne  percentage  of  each  constituent  in  the 
mixture. 

Table  2 

DAMAGE  THRESHOLDS  FOR  VAPOR- PI iASE  MIXTURES 
(adapted  from  reference  14) 

N 9.05  x 103 

n2-l  /as(n2-l ) 


Zns 

_ ThFu 

Index 

fra  d (MV/ cm) 

(1028  m"3) 

MV/cm 

0 

1QC 

1.5 

0.65 

4.92 

4.55 

23 

77 

1.59 

0.429 

3.86 

3.69 

3', 

65 

1.72 

0.406 

2.97 

2.88 

50 

50 

1.91 

0.34 

2.12 

2.12 

80 

20 

2.15 

0.223 

1.46 

1.53 

90 

10 

2.33 

0.223 

1.17 

! .25 

0 

0 

tr~  ‘ 

0 

2.36 

0.175 

1 .11 

1.21 

66 


Figure  17  Experimental  Threshold  Field  versus  Predicted  Threshold  Field 
of  Differing  Thickness 


AFWL-TR- 76-61 


These  data  are  plotted  in  figures  19  and  20.  In  each  case  the  fit  to  a least 
square  is  very  good,  being  ±7.5  percent  in  the  former  and  ±7.2  percent  in  the 
latter. 

The  next  previously  published  data  considered  are  those  collected  by 
Fradin  and  co-workers  (ref.  17).  They  damaged  the  bulk  of  11  alkali  halves 
at  1.06  ym.  Table  3 is  a compilation  of  the  data. 


Table  3 

BULK  DAMAGE  THRESHOLDS  FOR  ALKALI  HALIDES  WITH  4.7  NSEC  PULSES 
(adapted  from  references  15,  17) 


Material 

Index 

Field  (MV/cm) 

(1028  m"3) 

9.05  x 1C3  rm 

Nal 

1.74 

1.14 

1.46 

2.48 

NaBr 

1.62 

1.54 

2.31 

3.22 

NaCl 

1.53 

2.3 

3.32 

4.02 

NaF 

1.32 

3.77 

10.  G3 

7.99 

KI 

1.64 

G.62 

1.32 

2. 85 

KBr 

1.54 

0.87 

2.U4 

3.64 

KC1 

1.49 

1.31 

2.54 

4.19 

KF 

1.36 

2.99 

6.12 

5.51 

Rbl 

1.63 

0.92 

1.21 

2.85 

RbBr 

1.54 

1.27 

1.78 

3.55 

RbCl 

1.48 

1.54 

2.32 

4.18 

These  data  are  plotted  In  figures  21  and  22,  The  least  squares  fit  to  the  full 
expression  gives  a deviation  of  ±20.1  percent,  while  the  fit  to  N/n*-l  Is  good 
to  ^21 . 9 percent,  the  latter  being  close  to  the  quoted  experimental  precision  of 
±20  percent  (ref.  17). 

Before  examining  the  pulse  length  dependence  of  laser-induced  damage,  it  is 
appropriate  to  ask  whether  the  processes  at  30  psec  are  the  same  as  the  damage 
processes  at  40  nsec.  The  data  collected  by  Smith  (ref.  54)  answer  this  ques- 
tion. Smith  determined  the  threshold  electric  field  for  a variet;/  of  substances 
using  a 30-psec  Nd+3  In  YAG  laser  operating  at  1.06  ym.  The  variation  of  the 
laser-induced  damage  threshold  as  a function  cf  material  tested  followed  the 
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Figure  20.  Experimental  Threshold  Field  versus  N/n2-!  for  Vapor-Phase  Mixture  Films 


Figure  21.  Experimental  Bulk  Threshold  Field  versus  Predicted  Field  for  Alkali  Halides 


Figure  22.  Experlaental  Bulk  Threshold  Field  versus  N/n2-l  for  Alkali  Halides 
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functional  form  given  In  equation  (62).  This  Is  taken  as  an  Indication  that 
the  damage  process  does  not  change  In  its  essential  aspects  as  the  pulse  dura- 
tion Is  shortened  from  4C  nsec  to  30  psec.  The  essential  difference  Is  that 
the  threshold  fields  are  nearly  an  order  of  magnitude  higher  at  30  psec  than  at 
40  nsec.  This  variation  will  be  tested  against  the  behavior  predicted  by 
equation  (76). 


First  consider  the  material  variation  of  the  laser-induced  damage  threshold 
at  30  psec  pulse  duration.  Smith's  data  (ref.  54)  are  collected  in  tablo  4 and 
presented  In  figure  23. 

Table  4 

SULK  DAMAGE  THRESHOLDS  FOR  ALKALI  HALIDES  WITH  30  PSEC  PULSES 
(adapted  from  reference  54) 


MateM  &1  Index 


Field  (MV/cm) 


9.05  x 10s 
^s(n*-1) 


NaF 

1.32 

10.77 

7.99 

NaCl 

1.53 

7.34 

4.01 

NaBr 

1.62 

5.67 

3.22 

KF 

1.36 

8.34 

6.31 

rci 

1.49 

5.86 

4.19 

KBr 

1.54 

5.33 

3.64 

KI 

1.64 

5.87 

2.85 

L1F 

1.38 

12.24 

7.04 

Rbl 

1.63 

3.40 

2.85 

The  measured  field  fits  the  dependence  predicted  by  the  final  column  In  table  4 
to  ±17.4  percent.  Thus,  the  basic  damage  mechanism  does  not  change  for  pulse 
variations  between  30  psec  and  40  nsec. 


As  the  next  topic  In  the  review  of  previously  published  data  consider  the 
pulse  length  dependence  predicted  by  equation  (76).  These  data  collected  by 
Smith  (ref.  54)  and  Fradln  (ref.  17)  can  l compared  because  they  used  the  same 
equipment  to  perform  their  experiments.  The  data  are  presented  In  table  5 and 
figure  24.  Table  5 lists  the  threshold  daiMge  field  for  each  material  at  the 
specified  pulse  durations.  The  NaCl  data  are  presented  In  figure  24.  The  least 
square  fit  Indicates  a dependence  of  threshold  field  versus  pulse  length  as  Eth  ■ 
Const  tp"0,257  which  comparts  favorably  with  E^h  • Const  tp"0,21  predicted  by 
equation  (76).  The  data  fit  the  line  to  ±8.5  percent  which  Is  remarkably  good. 
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Table  5 

BULK  DAMAGE  THRESHOLDS  FOR  FIVE  ALKALI  HALIDES  AND  FUSED  SILICA 
(adapted  from  references  16,  17,  54) 


10.3  nsec 

4.7  nsec 

300  psec 

30  psec 

NaCI 

2.1  (MV/cm) 

2.3 

4.7 

7.34 

NaT 

3.77 

10.77 

NaBr 

1.54 

5.57 

KF 

2.99 

8.34 

Rbl 

0.92 

3.40 

Fused  Silica 

5.2 

11.68 

As  the  final  comparison  between  the  theory  and  previously  published  results 
consider  the  spot  size  dependence  reported  by  DeShazer  et  al.,  (ref.  55).  They 
damaged  the  surface  of  single  crystal  rutile  (T102).  The  focused  spot  size  was 
varied  from  50  urn  to  200  um  by  using  a selection  of  lenses  with  various  focal 
lengths.  The  data  are  presented  In  figure  25  as  a plot  of  threshold  Intensity 
versus  inverse  spot  size.  The  agreement  between  the  experimental  points  and  - 
the  relationship  predicted  by  equation  (78)  Is  11.1  percent  over  this  range  of 
spot  sizes. 

6.  RESULTS  OF  THE  CURRENT  STUDY 
a.  Thin  Films 

Surfaces  of  optical  materials  do  not  exhibit  damage  thresholds  as  high 
as  their  bulk  values.  To  the  author1 $ knowledge  there  has  been  only  one  report 
which  claimed  that  the  surface  and  bulk  threshold  of  a material  were  measured 
to  be  the  same  (ref.  16).  It  was  shown  (ref.  9)  that  clean  optical  surfaces 
exhibited  higher  damage  thresholds  than  do  dirty  surfaces.  In  fact,  a piece  of 
dust  on  an  optical  surface  can  absorb  laser  radiation,  and  the  resulting  tem- 
perature rise  Is  sufficient  to  cause  localized  damage.  It  was  also  shown 
(ref.  56)  that  particular  polishing  compounds  containing  material  which  Is 
highly  absorbent  at  the  laser  wavelength  will  lower  the  damage  threshold.  Spe- 
cifically, jeweler's  rouge  (iron  oxide)  lowers  the  surface  threshold  for 
1.06-um  laser  wavelength  even  though  the  Iron  oxide  is  not  present  In  sufficient 
quantity  to  be  Indicated  by  Auger  spectroscopy.  However,  even  If  these  two 
problems  of  dirt  and  particulate  Inclusions  are  scrupulously  avoided,  there  still 
remains  the  problem  of  surface  roughness.  A major  portion  of  an  -*ffort  parallel 
to  the  one  reported  here  (ref.  50)  was  concerned  with  the  evaluation  of  surface 
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Figure  25.  Threshold  Intensity  versus  Focused  Spot  Size  for  TiO 
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roughness  as  It  affects  surface  thresholds  to  laser-induced  damage.  In  that 
study  a series  of  fused  silica  sanples  with  root  mean  square  surface  roughness 
varying  over  nearly  two  decades  was  prepared  In  triplicate.  One  of  each  rough- 
ness was  lef*  uncoated,  one  of  each  roughness  was  coated  with  a half  wavelength 
thickness  cf  S1QZ,  and  one  of  each  roughness  was  coated  with  a half  wavelength 
tnlckness  of  MgFz.  The  results  of  that  study  Indicated  a damage  behavior  given 
by 

om  E^  ■ constant  (82) 

for  the  three  experiments.  The  value  of  m was  0.61  for  the  uncoated  samples, 
0.455  for  the  MgFz  coated  samples,  and  0.42  for  the  S102  coated  samples.  So 
that  the  Intrinsic  strength  of  several  optical  surfaces  can  be  compared.  It  Is 
necessary  to  eliminate  the  confusing  effects  of  varying  surface  roughnesses. 

This  can  be  done  by  using  equation  (82)  to  compare  the  threshold  breakdown 
fields  at  a standard  roughness.  To  do  so,  let  be  a standard  si  rface  rough- 
ness for  a particular  set  of  samples.  Then  If  o^  Is  the  surface  roughness  of 
the  1th  sample  In  the  set,  the  threshold  field  for  sample  1,  If  it  were  to  have 
the  standard  roughness  o^.  Is  given  by 

EJh  ’t«k>  • (£)"  E«<  (0i>  (83) 

In  that  which  follows,  m Is  taken  as  0.5,  the  average  of  the  experimentally 
determined  values.  For  the  thin  film  data  to  be  presented,  was  arbitrarily 
taken  as  9.8  K rms,  the  roughness  of  the  PrF,  coated  sample,  while  for  the  bare 
surfaces  a ak  of  13.5  ! rms  Is  used,  it  being  the  MgF2  roughness. 

The  thln-fllm  data  are  presented  In  table  6 and  figures  26  and  27.  Again 
the  threshold  field  Is  plotted  against  both  N/n*-l  and  9.05xl0*//as(n2-l ) . 

The  films  marked  with  an  asterisk  proved  Inhomogeneous  In  their  spectral 
reflectivity  curves.  As  such,  they  are  expected  to  behave  In  less  than  an  Ideal 
manner.  Figure  26  Is  a plot  of^c  Eth  versus  9.05xl0*//os(n2-l).  Tha  three 
Inhomogeneous  films  are  Indicated  by  plus  signs  to  mark  their  positions.  All 
three  fell  below  a least  squares  fit  made  without  their  Inclusion.  The  perfor- 
mance of  LI F and  MgFt  in  the  bare  surface  data  to  follow  will  be  particularly 
revelatory  In  line  with  the  contention  that  both  LI F and  MgF2  should  exhibit 
considerably  higher  damege  thresholds  than  did  their  half-wave  films.  A similar 
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Table  6 

ROUGHNESS  CORRECTED  THRESHOLDS  FOR  10  HALF-WAVE  FILMS  ON  FOUR  SUBSTRATES 

N 9.05xl03 


Material 

Index 

a(A) 

a° • s Fi el d 
(MV/cm) 

n2-l 

(102V3) 

/as(n2-l 

(MV/cm) 

LiF*  on  F3 

1.38 

12.34 

0.578 

13.356 

7,04 

MgF2*  on  F5 

1.377 

17.08 

1.15 

9.705 

6.72 

Si02  on  FS 

1.449 

17.4 

1.18 

6.02  . 

5.24 

S1O2  on  BK-7 

1.449 

20.8 

1.09 

6.02 

5.24 

3aF2*  on  F5 

1 .45 

12.3 

0.554 

4.515 

4.98 

Ti.F  on  FS 

1.49 

12.3 

0.990 

5.045 

4.66 

CeT3  on  FS 

1.55 

12.9 

0.71 

4.636 

4.09 

Prt-3  on  FS 

1.62 

9.8 

0.906 

6.216 

3.80 

A1203  on  Sapphire 

1.754 

15.9 

0.556 

5.58 

3.00 

Zr02  on  FS 

2.0 

15.1 

0.48 

2.688 

1.98 

ZnSe  on  ZnSe 

2.40 

31.5 

0.207 

0.9S 

1.13 

study  of  BaF2  was  not  possible  due  to  the  lack  of  a bulk  sample  of  BaF?.  The 
fit  is  very  good  for  the  points  considered  (i.e.,  all  but  the  t ee  inhomogene- 
ous films),  having  an  average  error  of  ±12.2  percent. 

Figure  27  is  a plot  of  the  roughness-corrected  threshold  field  versus 
N/n2-l.  Here  the  fit  is  not  quite  so  good  at  ±19.4  percent;  however,  the  gen- 
eral trend  is  still  apparent. 

From  both  figures  an  experimental  justification  can  be  seen  for  the  use 
of  the  -oughness  to  obtain  "equivalent"  surfaces.  From  the  raw  data  the  uncor- 
rected threshold  fields  for  the  Si 02  films  on  fused  silica  and  BK-7  were  0.892 
and  0.757  MV/cm,  respectively,  on  surfaces  with  17.4  and  20.8  A rms  roughnesses. 
The  correction  for  surface  roughness  reconciled  the  two  figures  to  within  8 per- 
cent which  is  within  the  experimental  accuracy. 

The  fact  that  a material  is  ir  film  form  usually  relegates  its  thresholc 
to  a value  below  that  of  a bare  surface  cf  the  same  material.  Films  suffer  fron 
problems  in  adhesion  different  mechanical  and  crystallographic  properties  from 
the  bulk  material,  fusions,  residua-  stress,  and  varying  degrees  of  Inhomoge- 
neity. Table  7 is  a comnilation  of  the  ratio  of  thin-rilm  to  bare  surface 


♦Denotes  Inhomogeneous  films,  FS  « fused  silica 
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thresholds  observed  during  this  investigation.  Both  the  thin-film  and  surface 

O 

thresholds  have  been  corrected  for  roughness  to  an  equivalent  13.5  A rms 
surface. 

Table  7 

RATIO  OF  THIN-FILM  TO  BARE  SURFACE  THRESHOLD 


Material 

Bare  surface  threshold 

Film  threshold 

Ratio 

Si02  on  FS 

2.09 

0.96 

0.459 

LiF  on  FS 

3.09 

0.49 

0.159 

MgF2  on  FS 

2.78 

0 98 

0.159 

A120_.  on  Sapphire 

1.09 

0.47 

0.43 

ZnSe  on  ZnSe 

0.142 

0.176 

1.26 

The  inhomogeneous  films  on  Li F and  MgF2  exhibit  the  lowest  ratios.  The  ZnSe 
film  is  interesting  in  that  it  has  a higher  threshold  than  the  bare  ZnSe  sur- 
face. This  result  if  not  completely  unexpected  due  to  the  difficulty  in 
producing  clear,  unbanded,  homogeneous  bulk  ZnSe  (ref.  57).  Even  though  one's 
first  thought  would  be  to  never  place  a dielectric  coating  on  a material  since 
films  exhibit  lower  thresholds  than  uncoated  surfaces,  there  are  good  reasons 
for  doing  so.  An  anti  reflective  coating  can  consist  of  a single  material  a 
quarter  wavelength  thick,  with  Index  na,  placed  on  a surface  with  refractive 

a 

index  nfa  in  which  na  *-^n^.  Thus  S102  with  index  1.449  would  make  a quarter- 
wave  anti  reflective  coating  for  a material  with  index  2.1.  A typical  threshold 
for  a material  with  n * 2.1  is  about  0.7  MV/cm  for  a 12.5  & rms  surface  finish. 

Thus,  S10  with  an  experimental  thin-fllm  threshold  of  0.96  MV/cm’ could  be 

2 

used  successfully  an  an  anti  reflective  coating  and  may  in  fact  increase  the 
damage  threshold  of  the  system  since  the  maximum  field  will  be  at  the  air-film 
Interface.  The  same  comments  hold,  of  course,  for  any  two  materials  in  the 
proper  ratio,  although  many  anti  reflective  coatings  are  placed  on  materials  with 
such  low  Index  that  a suitable  material  Is  not  available  with  the  proper  index 
to  make  a quarter-wave  anti  reflective  coating.  In  such  cases  multilayer  anti- 
reflective  coatings  are  used.  To  show  that  the  weakest  layer  in  a multilayer 
coating  will  fall  when  the  electric  field  reaches  Eth,  a series  of  multilayer 
films  were  subjected  to  laser  damage. 
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Two  types  of  film  stacks  were  tested:  alternating  even-numbered  half- 
wave layers  of  Zr02  and  Si02  and  alternating  odd-numbered  half-wave  layers  of 

ZrO  and  SiO  . The  even-numbered  series  consisted  of  2,  6,  and  10  layers  with 

2 2 

the  Zr02  outermost.  The  odd-numbered  series  consisted  of  3,  7,  11,  and  15 
layers;  each  having  Zr02  at  the  inner  (i.e.,  on  the  fused  silica  substrate)  and 
outer  layers.  Table  8 is  a compilation  of  the  data  gathered  on  these  film 
stacks. 


Design 

THRESHOLD 
Threshold  Energy 

(mj) 

Table  8 

OF  FILM  STACKS 

Incident  Field 
(MV/cm) 

Threshold  Field 
(MV/  on) 

g(LH) 

9.7 

1.03 

0.53 

g (LH ) 3 

4.53 

0.70 

0.62 

g(LH)s 

3.6 

0.63 

0.60 

gH(LH) 

3.17 

0.59 

0.48 

gH (LH ) 

3.5 

0.62 

0.58 

gH(LH)5 

3.78 

0.64 

0.61 

gH(LH) 7 

4.22 

0.68 

0.66 

Here  g denotes  the  substrate,  H the  high  index  film  (Zr02)  and  L the  low  index 
film  (Si02).  In  support  of  the  field  viewpoint  of  laser-induced  damage  it  is 
the  field  In  the  Zr02  which  is  the  most  nearly  constant.  The  threshold  field 
is  0.58  MV/cm  ±10.4  percent  while  the  threshold  energy  is  4.64  mj  ±49  percent 
and  the  incident  optical  field  Is  0.699  MV/cm  ±21.6  percent. 

In  an  attempt  to  discover  the  best  method  of  film  deposition,  MgF2 
films  were  placed  on  fused  silica  substrates  by  radio  frequency  sputtering, 
electron  beam  heating,  and  thermal  evaporation  techniques.  The  resultant 
threshold  fields  are  given  in  table  9. 

Table  9 

THRESHOLDS  FOR  MgF2  FILMS  VERSUS  DEPOSITION  TECHNIQUE 

Technique  Threshold  Electric  Field 

Electron  Gun  0.737  MV/cm 

Thermal  Evaporation  0.797  MV/cm 

R-F  Sputter  0.899  MV/cm 
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At  best  the  results  are  inconclusive.  The  variables  involved,  both  those  known 
and  thos  : unsuspected,  are  not  we’l  understood.  The  technique  which  works  well 
for  one  material  will  not  necessarily  produce  a good  film  with  another  material. 

The  final  result  to  be  discussed  in  thin  films  was  due  to  a bit  of 
serendipity.  A set  of  thin-film  samples  of  various  materials  had  been  prepared 
by  electron  beam  heating  with  the  angle  of  Incidence  of  the  coating  material 
as  the  variable.  When  the  data  were  analyzed  they  made  no  sense.  True,  there 
were  wide  variations  in  surface  thresholds,  but  there  was  not  a clear  trend 
with  angle  of  Incidence.  The  surface  roughnesses  of  the  substrates  were  checked 
but  proved  Inadequate  to  resolve  the  discrepancies.  However,  while  checking 
the  absorption  spectra  which  are  routinely  run  on  each  coating  as  a quality 
control  measure,  an  absorption  peak  at  1.405  pm  was  found  to  vary  in  direct  or 
inverse  proportions  with  the  damage  threshold,  depending  on  the  film  material. 

A check  was  made  which  indicated  that  the  fused  silica  substrate  had  not  all 
been  supplied  to  the  coating  company  by  the  same  vendor.  In  fact,  the  absorp- 
tion peak  was  found  to  correspond  to  OH"  absorption  and  was  in  direct  proportion 
to  the  values  of  OH"  concentration  quoted  in  the  vendors'  literature  (ref.  29). 
The  data  presented  in  taole  10  summarize  the  experimental  results.  The  results 
are  quite  Interesting,  especially  in  the  S102  films.  In  that  material,  thresh- 
old field  varies  directly  with  OH"  concentration  In  the  substrate.  This  may  be 
because  Si 0^  goes  onto  a substrate  not  fully  oxidized  in  S10  form.  The  addi- 
tional oxidation  by  the  OH"  of  the  substrate  improves  the  film  (ref.  29).  Con- 
versely, in  the  MgF2  films  the  additional  oxidation  by  higher  OH"  concentration 
(seen  by  the  lower  transmission  at  1.4  pm  and  2.23  pm)  causes  the  MgF2  to  de- 
posit at  least  near  the  surface  as  MgO  interspersed  with  MgF2  (ref.  29).  The 
trend  is  suggested  in  A1203  films  for  improved  films  by  additional  oxidation  and 
no  clear  relationship  Is  evidenced  by  the  Zr02  films.  A definitive  experiment 
with  a controlled  set  of  samples  is  obviously  called  for  by  these  preliminary 
results. 

b.  Bare  Surfaces 

A set  of  11  uncoated  samples  was  tested.  The  samples  ranged  from 
crystalline  materials  such  as  L1F  and  sapphire  to  optical  glasses  such  as  BK-7 
and  LASr-6.  A wide  range  of  refractive  index  was  of  prime  concern,  and  the 
materials  tested  met  this  goal  with  Indexes  from  1.377  to  2.485.  Surface 

O 

roughness  was  specified  as  best  finish  and  ranged  from  13.5  A rms  for  MgF2  to 
80  for  L1F.  The  latter  was  a loan  from  a group  Interested  in  testing  bulk 
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Table  10 


DAMAGE  THRESHOLDS  FOR  X/2  FILMS  VERSUS  OH"  CONCENTRATION  IN  SUBSTRATE 


Material 

Threshold  field 
(MV/cm) 

* T @ 1 . 4 pm 
(%) 

Angle  of 
vapor  incidence 
(deg) 

Al2°3 

0.408 

90 

20.3 

ai20, 

0.391 

81 

29.3 

ai2o3 

0.315 

90 

42.0 

SiO 

2 

1.24 

83 

24.96 

CO 

o 

N>  1 

1.16 

85 

20.3 

S10 

2 

0.813 

92.5 

42 

S102 

0.681 

93 

29.3 

SiO, 

2 

0.564 

94 

18.55 

MgF2 

0.61 

95 

20.3 

MgF2 

0.50 

85 

18.55 

MgF2 

0.50 

85 

24.96 

MgF2 

0.50 

86 

29.3 

MgF2 

0.50 

86 

42.0 

Zr02 

0.198 

79 

20.3 

Zr02 

0.252 

80 

42.0 

ZrO, 

0.251 

85 

24.96 

ZrO 

2 

0.197 

86 

29.3 

ZrO 

2 

0.292 

87 

18.55 

properties  of  Li F and  as  such  no  emphasis  was  placed  on  obtaining  a highly 
polished  surface,  although  50  X rms  Is  about  the  best  that  can  be  achieved  on 
L1F  (ref.  58).  The  surface  threshold  field  has  been  corrected  for  roughness 
so  that  each  value  Is  equivalent  to  13.5  X rms  surface.  The  data  are  presented 
In  table  11  and  figures  28  and  29.  These  data  are  plotted  in  figures  28  and 
29.  The  agreement  with  equation  (62)  is  generally  good  with  the  exception  of 
two  of  the  glasses  in  figure  28.  The  calculated  fit  of  ±19.3  percent  is  near 
the  experimental  accuracy.  The  performance  of  SF-12  and  SF-14  is  unexplained. 
The  fit  without  these  points  becomes  ±9.5  percent.  The  parameters  for  these  two 
glasses  were  calculated  from  data  supplied  by  the  vendor*.  The  roughness 


*0r.  Karl  H.  Mader,  Schott  Optical  Company. 
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Table  11 

ROUGHNESS  CORRECTED  SURFACE  DAMAGE  THRESHOLDS  FOR  11  DIELECTRICS 


N 9.05x10* 

n2'-l  /05  (n*-l ) 

Material  Index  ga,s  Field  (MV/cm)  (1 0a8m~ 3 ) (MV/cm) 


L1F 

1.38 

3.09 

13.356 

7.04 

MgF2 

1.377 

2.78 

9.705 

6.72 

Si02 

1.449 

2.09 

6.02 

5.24 

BK-7 

1.507 

1.80 

5.96 

4.63 

A1  0 . 
2 j 

1.754 

1.09 

5.12 

3.0Q 

SF-12 

1.63 

0.75 

3.99 

3.48 

SF-14 

1.734 

0.61 

3.98 

2.88 

LASF-6 

1.933 

0.80 

5.01 

2.19 

—1 

O 

N 

2.28 

0.347 

2.29 

1.46 

ZnS 

2.288 

0.213 

1.194 

1.30 

ZnSe 

2.485 

* 0.142 

0.874 

1.04 

corrected  threshold  electric 

field  Is  plotted 

against  N/n2-1 

in  figure  : 

the  agreement  suffers  being  only  ±23.3  percent. 

A final  bit  of  interesting  Information  can  be  extracted  *rom  the  study 
by  House  (ref.  50).  Table  12  compiles  these  data. 

Table  12 

BREAKDOWN  FIELD  VERSUS  DISORDER 


Material 

Disorder  (0) 

Threshold  field 

9.05xl03 

/cr(n2-l) 

sio2 

(Bulk) 

2.53  % 

5.2  MvW163 

5.24  MV/cm 

S10 

(Surface) 

19.8  (FECO) 

1.73 

1,85 

</> 

O 

PO 

(Surface) 

42.0  (FECO) 

0.987 

1.27 

SiO, 

z 

(Surface) 

1*0.0  (TIS) 

0.359 

0.695 

SIO 

2 

(Surface) 

773.4  (TIS) 

0.322 

0.296 

S10s 

(Surface) 

926.5  (TIS) 

0.240 

0.270 

This  Indicates  a definite  relationship  between  disorder  In  the  material  and  the 
threshold  field  for  laser-1 nduced  breakdown.  For  lack  of  a better  term,  dis- 
order was  chosen  and  taken  as  the  average  atomic  spacing  In  the  bulk,  o , 
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Figure  29.  Roughness  Corrected  Surface  Threshold  Field  versus 
N/n*-l  for  11  Dielectric  Materials 
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and  the  root  mean  square  surface  roughness,  o,  for  those  samples  damaged  on 
the  surface.  The  values  of  surface  roughness  for  o less  than  100  A rms  were 

O 

obtained  by  FECO  Interferometry,  while  those  over  100  A rms  were  obtained  by 
the  method  of  total  Integrated  scatter.  FECO  Is  not  effective  above  100  A rms 
according  to  Dr.  J.  Bennett  of  the  Michel  son  Laboratory  at  the  Naval  Weapons 
Center  of  China  Lake,  California  (ref.  59).  The  fields  quoted  are  not  cor* 
rected  for  roughness  but  are  the  actual  measured  values.  The  threshold  fields 
are  compared  with  calculations  from  equation  (62)  in  the  final  column,  using 
the  disorder  in  place  of  c$  which  appears  In  equation  (62).  The  data  fit  the 
predicted  dependence  to  ±21.0  percent  over  two  dec*  '^s  of  the  variable.  The 
intriguing  point,  of  course,  is  that  the  fit  of  th*  data  is  not  absurdly  bad. 
Additional  calculations  are  presented  In  table  13  In  which  experimentally 
determined  damage  fields  are  used  in  equation  (62)  to  predict  disorders  which 
are  compared  with  measured  (surface)  or  calculated  (bulk)  values. 


Table  13 

CALCULATED  VERSUS  MEASURED  DISORDER 


Material 

Threshold 

field 

Disorder 

9.05xl03 

ToTn^TT 

Calculated 

disorder 

MgFz  (Surface) 

2.78  MV/cm 

13.5  A 

2.75  MV/cm 

13.2  A 

L1F  (Surface) 

1.27 

80 

1.19 

62.1 

Sapphire  (Surface) 

0.608 

43.0 

0.66 

51.4 

BK-7  (Surface) 

1.446 

21.0 

1.55 

24.2 

BK-7  (Bulk)  (ref.  16) 

4.7 

2.36 

4.63 

2.29 

LASF-6  (Surface) 

0.8 

19.94 

0.74 

17.1 

T10j  (Surface) 

0.285 

20.0 

0.48 

57.2 

ZnS  (Surface) 

0.107 

53.56 

0.29 

399 

ZnSe  (Surface) 

0.116 

20.15 

0.39 

227 

Here  again  the  disorder  Is  the  average  atomic  spacing  In  the  bulk  and  the  root 
mean  square  roughness  for  the  surface  studies.  The  six  cases  considered  In 
table  12  and  the  first  six  In  this  table  give  good  agreement  between  threshold 
field  and  the  value  predicted  by  using  equation  (62)  with  a fit  of  ±16.2  percent. 
The  calculated  dlsordtr  Is  obtained  by  Inserting  the  measured  value  of  Eth, 
threshold  field.  Into 


. 9.05  x 10s 
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SECTION  V 
CONCLUSIONS 

The  conclusions  for  the  bare,  uncoated  optical  components  Is  fairly 
straight-forward.  To  maximize  the  power  density  which  can  be  transmitted 
through  an  optical  component,  one  needs  to  maximize  the  quantl tyj-^o^  (n*-l)J‘l. 
Other  parameters  are,  of  course  Important  but  to  a lesser  degree.  For  example, 
linear  absorption  Is  generally  higher  In  high  Index  materials  (ref.  31).  This 
contention  would  go  a long  way  toward  explaining  why  the  three  high  Index 
materials  (T102,  2.28;  ZnS,  2.288;  ZnSe,  2.485)  showed  such  poor  performance 
In  relation  to  the  threshold  fields  predicted  by  equation  (62).  If  laser- 
induced  damage  Is,  In  the  end,  a process  of  thermal  distortion,  melting,  or 
crazing,  then  any  mechanism  which  more  effectively  couples  the  laser  energy 
Into  the  lattice  affords  a means  for  lowering  the  damage  threshold.  One  must 
then  conclude  that  an  increase  In  linear  absorption  will  result  In  a decrease 
In  damage  threshold.  This  contention  is  borne  out  by  evidence  collected  during 
the  course  of  this  study.  Several  samples  of  two  glass  types  were  supplied  with 
as  nearly  equal  surface  finish  as  possible.  The  glass  types  were  Owens-Illinois 
ED-4  and  ED-2.  The  only  difference  in  the  two  glasses  is  that  ED-4  Is  undoped 
and  ED-2  is  a Nd+3  doped  laser  glass.  The  refractive  indexes  are  both  the  same 
and  the  number  densities  are  essentially  the  same.  The  physical  difference 
is  that  the  doped  ED-2  exhibits  an  Increased  linear  absorption  at  the  laser 
wavelength  of  1.06  urn.  The  matched  pairs  of  ED-2  and  ED-4  were  finished  in  a 
variety  of  ways  utilizing  different  conventional  polishing  and  Ion  polishing 
techniques.  In  each  matched  set  the  doped  ED-2,  with  a linear  absorption 
coefficient  of  only  0.001  to  0.002  cm"1  (ref.  60),  damaged  at  a lower  threshold 
than  did  the  undoped  ED-4  which  had  lower  absorption.  Facilities  were  not 
available  during  the  course  of  this  study  to  measure  absorption  coefficients  as 
low  as  0.002  cm'1.  The  three  high  Index  materials  in  table  13,  however,  were 
measured  for  transmission  using  a Carey  14  Spectrophotometer  and  in  two  of  the 
substances  high  absorption  was  found.  The  T102  had  a transmission  of  71.2  per- 
cent compared  to  the  71.9  percent  expected  from  a nonabsorbing  sample  with  an 
Index  of  2.28.  The  ZnS  had  a transmission  46.5  percent  compared  to  the  expected 
value  of  71c 7 percent  and  the  ZnSe  had  50.5  percent  compared  to  the  67  percent 
for  a nonabsorbing  material  with  an  Index  of  2.485.  From  table  13  one  can 
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compare  the  r?  os  of  the  measured  threshold  fields  to  the  field  predicted  by 
equation  (62).  For  the  three  materials  they  are:  Ii02,  0.59;  ZnS,  0.37; 

ZnSe,  0.30.  The  apparent  absorption  losses  for  the  three  are  0.7  percent, 

25.2  percent  and  16.5  percent,  respectively,  although  the  accuracy  of  the 
determination  does  not  allow  for  much  credence  in  the  0.7  percent  value.  Thus, 
the  apparent  trend  is  that  increasing  absorption  results  in  a lower  ratio  of 
the  observed  threshold  compared  to  that  predicted  by  equation  (62).  That  ZnSe 
does  not  perform  optimally  can  be  Inferred  from  the  fact  that  the  ZnSe  film 
exhibited  a higher  threshold  than  the  uncoated  ZnSe  surface.  Although  the  study 
of  laser-induced  damage  is  not  a sufficiently  detailed  science  to  allow  direct 
comparison  of  results  between  two  separate  Investigations,  the  work  of 
DeShazer  et  a!.,  (ref.  11)  suggests  strongly  that  both  T102  and  ZnS  may  be 
stronger  in  quarter  wave-film  than  the  uncoated  materials.  Intrinsic  absorption 
Is  not  the  only  mechanism  which  would  explain  the  anomalously  low  results  for 
the  high  Index  materials.  High  index  materials  tend  to  be  mechanically  soft. 
They  are  generally  hard  to  polish  in  that  the  polishing  compound  scratches  and 
pits  the  surface.  The  polishing  compound  can  also  easily  become  Imbedded  In 
the  material.  If  a relatively  large  amount  of  polishing  compound  becomes  im- 
bedded in  the  surface,  a simple  heat  absorption  mechanism  at  the  particulate 
inclusions  could  lower  the  threshold  field. 

Aside  from  the  high  index  materials  the  equation 

E a 9.05  x 103 

th  >£T  (n2  - 1)  (85) 

adequately  describes  the  observed  electric  field  threshold  for  laser-induced 
damage  to  various  optical  surfaces.  It  also  appears  to  be  predictive  of  the 
electric  field  threshold  for  laser-induced  bulk  damage  to  transparent  dielec- 
trics. Before  this  statement  can  be  made  conclusive,  a large  number  of  other 
materials  needs  to  be  tested  against  the  predicted  behavior. 

Finally,  consider  the  threshold  fields  for  the  bulk  of  dielectric  materials. 
From  figures  21  and  23,  It  is  apparent  that  the  relative  values  for  bulk  thresh- 
old fields  follows  the  dependence  suggested  by  equation  (62).  It  Is  also  clear 
that  the  damage  processes  are  essentially  the  same  at  30  psec  and  10  nsec.  The 
difference  In  the  absolute  values  of  the  field  between  15  psec  and  10  nsec 
can  be  explained  by  considerations  of  thermal  heating  of  the  plasma  \nd  lattice. 
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Thus,  the  proposed  model  explains  the  essential  variations  of  threshold  field 
as  a function  of  material  and  surface  properties  (both  on  the  surface  and  in 
the  bulk)  and  as  a function  of  laser  pulse  length. 

The  general  conclusion  for  thin  films  is  the  same  as  that  for  bare  surfaces. 
As  long  as  a thin  film  accurately  models  the  bulk  material,  equation  (62)  can 
be  used  to  predict  relative  electric  field  thresholds  between  different  film 
materials.  The  major  problem  involved  In  attempting  to  predict  thin-film 
thresholds  is  the  large  number  of  variables  Involved  in  the  coating  process. 
Residual  gas  pressures  in  coating  chambers,  temperature  of  the  substrate, 
deposition  rate,  deposition  technique,  and  condition  of  the  substrate  are  all 
important  considerations  In  thin-film  coatings.  The  deposition  of  thin  films 
appears  to  be  as  much  an  art  as  it  is  a science. 

The  major  specific  conclusion  about  thin  films  Is  that  the  roughness  of 
the  substrate  is  just  as  Important  in  determining  chin-film  thresholds  as  bare 
surface  thresholds.  No  longer  can  surface  or  thin-film  thresholds  be  given 
proper  credence  unless  the  surface  roughness  is  also  stated. 

Specific  material  properties,  such  as  free  radical  (OH-)  concentration, 
has  been  shown  to  be  important  In  some  thin-film  systems.  A factor  of  Z.Z  in 
electric-field  threshold  obtained  by  apparently  varying  only  the  type  of  fused 
silica  used  as  substrate  is  most  important  in  thin-film  technology. 

The  determining  factor  In  lase*"- induced  damage  is  the  optical  electric  field. 
Multilayer  thin-film  systems  can  be  increased  in  damage  resistance  by  utilizing 
this  fact.  Optical  systems  can  be  designed  to  minimize  electric-field  amplitudes 
at  ^eak  points  in  the  systems. 

The  proposed  model  does  not  deviate  strongly  from  the  observed  behavior. 
Electron  avalanche  was  not  considered  nor  iis  It  apparently  necessary  to  describe 
laser-induced  damage  of  optical  components.  Direct  local-field  generation  of 
free  electrons  can  account  for  the  observed  dependence  of  threshold  fields  on 
refractive  index,  number  density,  and  surface  roughness.  It  Is  seen  th«t  no 
separate  mechanism  need  be  postulated  for  surface  versus  bulk  damage  since  the 
threshold  fields  for  both  are  shown  to  obey  the  same  descriptive  relationship. 

In  conclusion,  a simple  model  cf  laser-induced  damage  has  led  to  a predictive 
formula  for  threshold  electric  field  which  holds  for  a wide  variety  of  materials 
In  bulk,  bare  surface,  and  thin-film  form. 
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